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Abstract 

We study the statistical properties of wave scattering in a disordered waveguide. The statistical 
properties of a "building block" of length 5L are derived from a potential model and used to 
find the evolution with length of the expectation value of physical quantities. In the potential 
model the scattering units consist of thin potential slices, idealized as delta slices, perpendicular 
to the longitudinal direction of the waveguide; the variation of the potential in the transverse 
direction may be arbitrary. The sets of parameters defining a given slice are taken to be statistically 
independent from those of any other slice and identically distributed. In the dense-weak-scattering 
limit, in which the potential slices are very weak and their linear density is very large, so that the 
resulting mean free paths are fixed, the corresponding statistical properties of the full waveguide 
depend only on the mean free paths and on no other property of the slice distribution. The 
universality that arises demonstrates the existence of a generalized central-limit theorem. 

Our final result is a diffusion equation in the space of transfer matrices of our system, which 
describes the evolution with the length L of the disordered waveguide of the transport properties 
of interest. In contrast to earlier publications, in the present analysis the energy of the incident 
particle is fully taken into account. 

For one propagating mode, N = 1, we have been able to solve the diffusion equation for a number 
of particular observables, and the solution is in excellent agreement with the results of microscopic 
calculations. In general, we have not succeeded in finding a solution of the diffusion equation. 
We have thus developed a numerical simulation, to be called "random walk in the transfer matrix 
space", in which the universal statistical properties of a "building block" are first implemented 
numerically, and then the various building blocks are combined to find the statistical properties 
of the full waveguide. The reported results thus obtained (in which use was made of a "short- 
wavelength approximation" ) are in very good agreement with those arising from truly microscopic 
calculations, for both bulk and surface disorder. 

PACS numbers: 05.60.Gg, 73.23.-b, 05.40.-a, 84.40.Az 
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I. INTRODUCTION 



The statistical theory of certain complex wave interference phenomena, like the statistical 
fluctuations of transmission and reflection of waves, is of considerable interest in many 
fields of physics [1-11]. In the literature one has contemplated situations in which such a 
complexity derives from the chaotic nature of the underlying classical dynamics, as in the 
case of chaotic microwave cavities and quantum dots, or from the quenched randomness of 
scattering potentials, as in the case of disordered conductors or, more in general, disordered 
waveguides. It is the latter domain that will interest us here. 

In studies performed in such systems one has found remarkable statistical regularities, 
in the sense that the probability distribution for various macroscopic quantities involves a 
rather small number of relevant physical parameters, or scaling parameters, while the rest of 
the microscopic details serves as mere "scaffolding" . In Ref. [12] it was shown that a limiting 
distribution of physical quantities indeed arises in the so called dense-weak- scattering limit 
(DWSL) and within a particular class of models: the individual, microscopic, scattering 
units were defined through their transfer matrices and an "isotropic" distribution of their 
phases was assumed. The limiting distribution that was found constitutes a generalized 
central-limit theorem (CLT). Within this model only one relevant physical parameter occurs: 
the mean free path, which is the only property arising from the individual scattering units 
that survives in the DWSL. This is consistent with the scaling hypothesis proposed by E. 
Abrahams et al. [13]. (When abandoning the DWSL, two parameters were needed in Ref. 
[14] to describe the conductance distribution.) The result found in Ref. [12] coincides with 
that of the maximum-entropy model that had been developed in Ref. [15], which gives 
rise to a diffusion equation known as the DMPK equation (after Dorokhov [16] and Mello, 
Pereyra and Kumar [15]), which can thus be interpreted as capturing the features arising 
from a CLT. 

CLT's associated with products of matrices had been studied earlier, as, for instance, in 
the well-known Oseledec theorem [6, 17]; the results of Refs. [12, 15] are consistent with 
this theorem in the localized regime. 

In spite of the successes of the diffusion equation of Ref. [15] in the study of the conduc- 
tance distribution [6], that equation fails to give the proper description when the difference 
in behavior of the various modes becomes relevant. A clear example was given in Ref. 
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[18], where the conductance distribution was studied in the crossover region (G) ~ e 2 /h. 
For waveguides with bulk disorder the description is excellent, whereas for waveguides with 
surface disorder it is not satisfactory [11, 19]. 

A class of limiting distributions wider than that of Ref. [12] was studied by one of the 
present authors and S. Tomsovic in Ref. [20] (to be referred to as MT), in which the isotropy 
assumption of Ref. [12] was relaxed to a large extent. The DWSL played an essential role 
and the result was a more general CLT than that of Ref. [12], expressed in terms of a 
generalized diffusion equation. The scaling parameters that appear in MT are the mean free 
paths (mfp's) for the various scattering processes that may occur in the problem. When the 
various mfp's can be represented by a single one, one encounters the same diffusion equation 
that was studied in Ref. [15] using a maximum-entropy model. Thus the MT model appears 
as a possible candidate to study, in the problem of waveguides with surface disorder, the 
influence of the specific scattering properties of the relevant modes. 

The ideas of MT are further developed in the "Brownian-motion" model of Ref. [21]: a 
waveguide of length L is enlarged by adding a piece of thickness SL [to be called a "building 
block" (BB)], small on a macroscopic scale but still containing many scatterers, which is 
likened to a Brownian particle which, in a time interval 5t, small on a macroscopic scale, 
still suffers many collisions from the molecules of the surrounding medium. The transfer 
matrix for a BB is written as M = I + e and the independent parameters which M depends 
upon are chosen so that e for the BB satisfies a number of properties, reminiscent of those 
of a Brownian particle: 

(e)sL = + O(5L 2 ) (1.1a) 
(bb)sl = 0(5L) , (1.1b) 

while higher moments of e behave as higher powers of 5L (see Eqs. (3.73) and (3.74) of Ref. 
[21]). The result of this analysis is the same diffusion equation as that of MT. 

Though appealing the assumptions behind the Brownian-motion model of Ref. [21] for 
the BB may be, they are, nevertheless, arbitrary. Of course, they can be deduced from the 
MT model for the more microscopic scattering units. However, even these have a certain 
degree of arbitrariness. It would be satisfactory if these models could be obtained in a 
unified way from a maximum-entropy "ansatz" : this, however, is not known to the present 
authors at this time. 
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The motivation of the present paper is to derive from a potential model the statistical 
properties of the BB and use them to find the "evolution" with length of the expectation 
value of physical quantities. Since the potential model will be introduced at the level of the 
individual scattering units, the approach to be presented here is, in a way, hybrid between 
the methods of MT and Ref. [21]. We shall see that within the present model it is not strictly 
true that the individual transfer matrices (resulting from the individual potentials) for the 
various scatterers are identically distributed [see Eq. (3.24) below], as was assumed in MT, 
and this fact will be taken into account. The continuous limit is also treated here in a more 
satisfactory way than in MT, and the energy appears explicitly in the following presentation, 
in contrast to earlier publications. We believe that the present model is physically more 
complete than that of Ref. [15]; it is also better founded than that of MT and Ref. [21], 
in the sense that there is a lesser degree of arbitrariness in the assumptions, although the 
resulting diffusion equation has a "structure" similar to the one obtained in MT and Ref. 
[21]. Our diffusion equation covers situations not contemplated in these references: we shall 
find it suitable to study wave-transport problems in which the physics of the various modes 
is relevant; we shall also find a good description of the statistical properties of quantities that 
involve phases, which were not described at all in previous models. The reader is referred 
to Ref. [22] for a preliminary account of the results of the present paper. 

The paper is organized as follows. In the next section we derive a Fokker-Planck equation 
for the "evolution" with the waveguide length L of the expectation value of the physical 
quantities of interest. That equation represents the central result of the present paper and 
is given in Eq. (2.10), which we reproduce here for convenience: 



We notice that Eq. (1.2) contains no drift term and is thus a diffusion equation. The physical 
observable is denoted by F(M), M being the transfer matrix of the sample of length L. The 
quantities D % ^ d (k, L) play the role of "diffusion coefficients", which are defined in terms of 
the second moments of e for the BB in Eq. (3.50) below (see the term linear in SL) and 
are given explicitly in Eq. (3.51) in terms of the mean free paths. Notice that the diffusion 
coefficients depend on the energy (~ k 2 ) and also on the length L of the sample. The mean 
free paths depend only on the second moments of the potential intensity of the individual 




ijhlXfi 
abcdafi 



ab,cd 




(1.2) 
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impurities [see Eq. (3.36)], higher moments being irrelevant for the diffusion equation: this is 
precisely what signals the existence of a CLT. In order to derive the diffusion equation (1.2) 
we need a statistical model for the building block (BB): this is derived in Sec. Ill using a 
potential model for the random impurities. Although the treatment of Sec. II is applicable 
to the orthogonal as well as to the unitary symmetry classes of Random-Matrix Theory 
(P = 1 and 2, respectively [23]), the potential model developed in Sec. Ill assumes time- 
reversal invariance, i.e., P — 1. Some of the specific relations derived there would have to be 
properly modified for the unitary case, P — 2. The results of Sec. Ill which are needed for 
the derivation of the diffusion equation (1.2) have an intrinsic interest as well, since they can 
be used to describe the statistical scattering properties of thin slabs. The diffusion equation 
(1.2) is first derived for arbitrary energy, and only later the short- wavelength approximation 
(SWLA) is contemplated; it is in this latter limit that some of the results obtained earlier 
can be recovered. Needless to say, we have no general way of finding either analytically or 
numerically the solutions of the above diffusion equation. We thus give in Sec. IV A some 
simple examples in which the analytic solution could be found; in Sec. IV B we develop a 
procedure to simulate numerically the diffusion process in transfer-matrix space, and present 
some of the results that we have been able to obtain so far. The conclusions of this work are 
given in Sec. V. A number of appendices have been included, in which some of the specific 
calculations are carried out. 



II. TRANSPORT IN Q-1D DISORDERED SYSTEMS. THE COMBINATION 
LAW AND THE SMOLUCHOWSKY EQUATION. THE DIFFUSION EQUATION 

Consider a q-lD disordered system of uniform cross section, connected, at both ends, 
to clean waveguides that support N open channels each. In the disordered region there 
is an underlying random potential to be specified later. In some applications we shall be 
concerned with a 2D waveguide with a width to be denoted by W. 

The scattering properties of the system will be described by means of its transfer matrix 
M, which can be written as 



M 11 


M 12 




a P 


M 21 


M 22 
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Each block M iJ (i = 1, 2; j = 1,2) in (2.1) is iV-dimensional, so that M is 2A r -dimensional. 
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(The block M 12 will occasionally be denoted by j3, a symbol not to be confused with the 
index for universality classes in Random-Matrix Theory.) One particular matrix element of 
the ij block will be designated as M^ b , where a, b (= 1, • • • , N) denote the channels. Some 
of the properties of the M matrix and its relation with the more conventional reflection and 
transmission amplitudes, which are elements of the S matrix, are summarized in App. A. 

The transfer matrix M will be considered to belong to one of the basic symmetry classes 
introduced by Dyson in Quantum Mechanics [23]. Here we shall be only concerned with 
scalar waves, so that, in applications to quantum mechanics, we shall only have "spinless 
electrons". In the "unitary" case, also denoted by (3 = 2, the only restriction on M is 
flux conservation, which is expressed by the pseudo-unit arity condition, Eq. (A4). In the 
"orthogonal" case (j3 = 1), time- reversal invariance imposes the restriction given by Eq. 
(A5). The "symplectic" case {(3 = 4) associated with half- integral spin will not be considered 
here. 

If the underlying potential has non-zero matrix elements between open and closed chan- 
nels, the 2iV-dimensional M matrix depends on an "effective potential" which contains 
information on closed channels, as explained in App. B. 

Consider now two non-overlapping scatterers. Their extended transfer matrices, Mi and 
M2 (which include open and closed channels, are infinite-dimensional and depend on the 
bare potential, as opposed to the effective one), have the multiplicativity property (see App. 
A3) 

M = M 2 Mi. (2.2) 

In past publications by one of the authors (PAM) (see, for instance, Ref. [10]), closed chan- 
nels have been neglected in the matrix multiplication of successive scatterers. In numerical 
simulations [11, 24] one sees that for individual configurations of the disordered system and 
in the calculation of the mean free path, the inclusion of closed channels is important. (In 
this paper, the expressions "closed channels" and "evanescent modes" will be taken as syn- 
onymous.) On the other hand, for the statistical fluctuations the conditions for neglecting 
the evanescent modes do not appear to be very stringent. For a given mean free path, the 
statistical properties of the different transport coefficients are found to be roughly indepen- 
dent of the number of evanescent modes (see also the discussion of the numerical simulations 
given in Sec. IV B 3). In this article we shall thus follow the earlier approximation and write 
the resulting transfer matrix as the product of the individual open-channel transfer matrices. 
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V 




FIG. 1: Schematic representation of a disordered wire and the building block (BB); (a) and (b) 
correspond to the different regimes (see Section III) defined by the inequalities given in Eqs. (3.1a) 
and (3.1b) (short-wavelength approximation (SWLA)), respectively. 

Suppose we start with a system containing n scattering units (to be denned at the begin- 
ning of Sec. Ill) and enlarge it by adding, on its right-hand side, say, a slab, to be called a 
building block (BB), containing m scattering units. Designating by the transfer matrix 
of the original system and by M( 5L *> that of the BB, the resulting transfer matrix is 

M (L+SL) = M (SL) M (L)^ ( 23 ) 

We assume the BB to be of arbitrary thickness 5L, and to contain many weak scatterers (see 
Fig. 1). Its transfer matrix M^ L "> will be written as 

M {5L) = I 2N + e. (2.4) 
The combination law for M, Eq. (2.3), can be written as 

M (L+6L) = M (L) + 6M (2.5a) 

= M {L) + eM {L \ (2.5b) 

Consider now a function F(M) of the M matrix, which we shall call an "observable": 
it could be, for instance, the transmission amplitude t a b, or the conductance G, which is 
proportional to the total transmission coefficient T. We are interested in the expectation 
value (F(M)) n of such an observable for a system containing n impurities. We first find 
below a recurrence relation with n for that expectation value and then, in the continuous 
limit, we shall find the equation that governs the "evolution" of (F(M)) L with increasing 
length L. 
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If we write a particular M matrix element as M^ h = £^ + irf^ b , we may consider the 
observable F(M) as a function of all those C a \ anc ^ rfab wri ich are relevant for the universality 
class (3 in question. For instance, in the orthogonal case ((3 = 1), because of the TRI relation, 
Eq. (A7), only the real and imaginary parts of the blocks M$ and M a \ 2 are relevant, while 
in the unitary case (j3 = 2) the real and imaginary parts of the four blocks M^ b (i = 1,2; 
j = 1,2) are needed. Alternatively, we may consider F(M) as a function of all the relevant 
M^'s and their complex conjugates, a possibility which will be found more convenient in 
what follows. For (3=1, Eq. (A7) shows that this is equivalent to expressing F(M) as a 
function of all the M^ b (i = 1,2; j = 1,2), enforcing Eq. (A7) at the end, while for (3 = 2 
we need all the M^ b (i = 1, 2; j = 1, 2) and their complex conjugates. In what follows we 
shall restrict ourselves to the orthogonal case. 

Writing the composition law for the M matrix as in Eq. (2.5), the expressions for 
the observable F(M) before and after adding the building block are related by the Taylor 
expansion 



F (M< n+m )) = F (AfW + SM) = F (M<">) + 



i\ UIV1 aa 
aa 



M=M( n ) 



aacf} 



+•••, (2.6) 



where the lower indices a, a, ■ ■ ■ on each M indicate channels and run over the values 
1, • • • ,N, while the upper indices i, A, • • • identify the block in Eq. (2.1) and take on the 
values 1,2. 

We take the expectation value of both sides of Eq. (2.6) with respect to the enlarged 
system containing n + m scattering units and use Eq. (2.5), considering the two pieces n 
and m to be statistically independent. We find 

(F(M)) n+m = (F(M)) n + £ <4> m (m£ 9 -^) 

ijX X M ' n 

aba 

+ iE('S«S.( J «^S;) +-■ < 2 - 7 > 

ijhlX)jb \ aa c(i I n 

abcda/3 

Here, (■ • •) denotes an average evaluated with the probability density for the transfer matrix 
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of the original sample containing n scattering units, i.e., 

(G(M))„e(G(mW)). (2.8) 

A few comments are in order at this point. We recall that the various matrix elements M^ b 
are not independent: they are related by the pseudounitarity condition, Eq. (A4), arising 
from flux conservation, and by Eq. (A5), associated with time-reversal invariance. Of course, 
we can express the M^'s in terms of independent parameters, like those occurring in the 
"polar representation" of the transfer matrix, and perform a Taylor expansion (similar to 
the one above) with respect to such independent parameters [10, 15]. Here, just as in Ref. 
[20] , we have found it simpler to perform the expansion with respect to the matrix elements 
M % l h themselves, as in Eq. (2.6), in the understanding that in the resulting expression (2.6) 
the M^'s and <5M*^'s have to be expressed in terms of independent parameters. The average 
appearing in Eq. (2.7) is thus performed with a probability density for such independent 
parameters. Since the latter depend on the underlying microscopic potentials, in this paper 
we shall not propose a model for the transfer matrix independent parameters, but rather for 
the microscopic potentials. (See also Ref. [20].) 

The next step is to describe the problem in the dense-weak- scattering limit (DWSL) briefly 
described in the Introduction (and defined in Eqs. (3.38) below), so that we can speak of 
the continuous length L of the system and the length 5L of the BB. Eq. (2.7) becomes 

(F(M)) L+SL = (F(M)) L + £ (4) /< d -^) 

ljX \ U1V1 aa I L 

aba 

+ S £ ^g^} + .... (2.9) 

abcda.fi 

To proceed, we need a statistical model for the BB. For this purpose, as we mentioned 
in the previous paragraph, a potential model is discussed in Sec. Ill, in which the BB is 
constructed as a collection of m individual scattering units represented by delta-potential 
slices. It is found that the first moment of e for the BB vanishes [see Eq. (3.29)], the second 
moments, in the DWSL, admit an expansion in powers of SL starting with SL itself (see Eq. 
(3.50), while higher moments behave as higher powers thereof (see the discussion following 
Eq. (D19)). Also, the very important result emerges that the dependence on the cumulants 
of the potential higher than the second drops out in the DWSL. These results are reminiscent 
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of the statistical behavior of the velocity increment of a Brownian particle during a time 
interval 5t during which many collisions from the surrounding medium have occurred [25]. 

When the moments of the BB, evaluated in the DWSL, are substituted in Eq. (2.9), we 
obtain, on the r.h.s of that equation, a power series in 5L. We also perform, on the l.h.s. 
of Eq. (2.9), a Taylor expansion of (F(M)) L+SL in powers of 5L around the "initial" value 
{F(M)) L . We can then identify the coefficients of the various powers of SL on the two sides 
of the equation. In particular, the coefficients of SL give the diffusion equation 



The quantities cd (k, L) play the role of "diffusion coefficients": they are defined in Eq. 
(3.50) below as proportional to the coefficient of the linear term in an expansion in powers 
of SL of the second moment of e for the BB and are given explicitly in Eq. (3.51) in terms 
of the mean free paths. The diffusion coefficients depend on the energy (~ k 2 ) and also on 
the length L of the sample. 

We remark that, just as the coefficients of SL in Eq. (2.9) are expressible in terms of 
the mfp's, the coefficients of higher-order terms in SL have a similar property, because the 
contribution of higher moments becomes irrelevant in the DWSL. Equating the coefficients 
of such higher-order terms on both sides of Eq. (2.9) we obtain results which could be 
derived from the diffusion equation (2.10) by successive differentiations. (See comment right 
after Eq. (3.55).) 

In the potential model discussed in the next section only the orthogonal case, f3 — 1, is 
contemplated. We expect a similar behavior for the unitary class, /3 — 2, although we do 
not have at the present moment the specific expression for each diffusion coefficient in this 
case. 

Eq. (2.10) represents the central result of the present paper. It depends only on the 
mean free paths which, in turn, depend only on the second moments of the individual delta- 
potential strengths [Eq. (3.36)]. The fact that cumulants of the potential higher than the 
second are irrelevant in the end signals the existence of a generalized CLT: once the mfp's 
are specified, the limiting equation (2.10) is universal, i.e., independent of other details of 
the microscopic statistics. 




(2.10) 



abcdaf3 
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III. STATISTICAL PROPERTIES OF THE BUILDING BLOCK 



In the present section we investigate the statistical scattering properties of the BB which 
was used in Sec. II to build a disordered system with a qlD geometry (see Fig. 1). 

Suppose that we model the scatterers constituting the BB by a sequence of thin slices 
(the scattering units referred to right above Eq. (2.3)) of cross section W D ~ X (D being; the 
dimensionality of the waveguide). From now on we denote the thickness of the slices by 2a 
and their separation by d. (See Fig. 2 below. Notice that in Fig. 1 the same symbols refer 
to individual scatterers; here, a slice may contain one or more of the individual scatterers 
shown in Fig. 1.) The statistical properties of the potential slices will be specified below (see 
Sec. IIIB1). Inside 2a, the r-th scattering slice is described by the potential V r (x,y). We 
denote by x the coordinate along the waveguide and by y the coordinates in the transverse 
direction. The distance d between slices is taken to be much larger than a, but much smaller 
than the wavelength A of the incident wave and the thickness 5L of the BB. Initially we 
do not specify the ratio of the wavelength A to 5L or the mean-free-path I (to be defined 
later), so we shall start out constructing the BB as a collection of m thin slices satisfying 
the inequalities 



Later on, in section HID, we shall find it advantageous to study a second regime, in which 
5L (and hence any final L) and I contain many wavelengths, i.e., 



corresponding to what we shall call the short-wavelength approximation (SWLA). 

In principle we have no restriction on the dimensionality D of the waveguide; however, 
to be specific, we shall restrict the discussion to two-dimensional waveguides with uniform 
width W. As we already indicated, in the potential model to be presented below we shall 
be concerned with the orthogonal, or — 1, symmetry class only. 

A. Properties of a single scattering slice 

Consider a single scattering slice with potential V(x,y) = h 2 U(x,y)/{2m), centered at 
the origin of coordinates x — 0, and let [C/(x)] a b be the matrix elements of U(x,y) with 



a < d < {\,5L,£}. 



(3.1a) 



a < d < A < {5L,£} , 



(3.1b) 
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respect to the "transverse" states Xa(y) of the waveguide, i.e., 

• w 



[U(x)] ab = I X a(y)U(x,y) Xb (y)dy, (3.2) 

with 



o 



x *m = \lw sin w> (3 - 3) 

a being an integer. Under the conditions 

ka < 1 (3.4a) 

K ab a < 1, (3.4b) 



where k = 2n / A = y2mE / h and 

K 2 ab = \U ab \ = \[U(0)] ab \, (3.5) 

we speak of a thin scatterer (a thin barrier or well) and the dependence of the potential 
across the thickness 2a is neglected. On the other hand, the quantity 

2a ■ U ab = u ah (3.6) 

(which has dimensions of k) is arbitrary. Such a scatterer can be well approximated by the 
"delta potential" 

U(x,y) = u(y)8{x), (3.7a) 
[U(x)} ab = u ab 5(x) , (3.7b) 

obtained formally taking the limits 

\U ab \ -> oo, (3.8a) 
a -> 0, (3.8b) 

in such a way that the quantity u ab of Eq. (3.6) stays fixed. From the inequalities (3.1) we 
see that the range 2a of the potential is the smallest length scale in the problem: the limit 
(3.8b) is the extreme idealization of this situation. 

Eqs. (3.7) define a delta-slice potential centered at the origin of coordinates. The poten- 
tial produced by the r-th delta slice, centered at X X ■y j IS written as 

U r (x,y) = u r (y)5(x - x r ) (3.9a) 

[Ur(x)] ab = {u r )ab S(X - X r ). (3.9b) 
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We remind the reader that U r (x,y) has dimensions of k 2 , whereas u r (y) and (u r ) a b have 
dimensions of k. 

A particle scattered by the potential of Eq. (3.7) inside the waveguide is described by 
the wave function 

oo 

^(x,y) = Y,Mx)lXa(y), (3.10) 

a=l 

which satisfies Schrodinger's equation; its components [f/H^L satisfy the coupled equations 
b + WWl = E MJ( X - *r) , 1 < a < iV (3.11a) 

^ ' 6=1 

( d 2 \ 00 

{Q^-rt)mx)l = E M 1 )]* ( u r) ab S(x - *r) , a>N + l. (3.11b) 

^ ' 6=1 

Eq. (3.11a) refers to open channels and Eq. (3.11b) to closed ones. The quantity k a , defined 
by the relation 



« = *"-(£)'. (3.12) 

is the "longitudinal" momentum for the open channel a, with the replacement fc a =>• m a for 
closed channels [10]. Notice that if Nit < kW < (N + the problem admits precisely N 
open channels. 

The open-channel, 2iV-dimensional, transfer matrix M (that relates open-channel am- 
plitudes on both sides of the potential) for the r-th slice, to be designated by M r , will be 
written as 



M r 



n M} 2 
121* rji/flir 



-?2Af + e r . (3.13) 



[Ml 2 ]* [Ml' 

Since, eventually, we shall be interested in the limit of weak scatterers in which M r is close to 
the unit matrix, we have introduced the difference e r between M r and the 27V-dimensional 
unit matrix I2N- In the above equation we have taken into account explicitly the fact that 
our system obeys time-reversal invariance (see App. A). The 11 and 12 blocks of the matrix 
e r are given by 

(e r )H = -i(v r Ue-^- k ^ ee (vrU(A)H (3.14a) 
feOS = -i{vrUe-^ +k ^ = (vr)aMH , (3.14b) 

where a and b label the open channels and thus run from 1 to N. We have defined 

= w*r))S = i(-ye^ jk ^-y +lk ^ (3.i5) 
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and we have introduced the real quantities 

(0 r U=5^=, (3.16) 

where, as explained in App. B, {u r ) a b is an "effective" potential strength that takes into 
account transitions to closed channels [see also Ref. [10], Eq. (3.134)]. 

In the above equations the strength of the various scatterers is arbitrary. As we already 
indicated, we shall be interested in the situation of weak scatterers, defined by the inequality 



(U r ab 



(3.17) 



which has to be added to the inequalities (3.1) in order to complete the specification of the 
physical regime. 

B. Construction of the Building Block. The regime (3.1a). 

1. The statistical model 

The BB is assumed, for the time being, centered at x = 0. For the application to Eq. 
(2.9) the BB will have to be translated to the interval (L, L + SL); this will be done in Sec. 
IIIC The BB is constructed from m delta slices located at the positions x r (see Fig. 2), 
i.e., assuming m to be odd, 

x r = rd (3.18a) 

m — 1 m — 1 , 

r = g — ,0,--. , — ^ — ( 3 - 18b ) 

SL = (m-l)d, (3.18c) 

where d denotes the distance between successive slices and SL the thickness of the BB. 

The m potentials u r {y), r = 1, ■ ■ ■ , m, are assumed to be statistically independent and 
identically distributed. We indicate the p-th moments of the individual M r (?/)'s and v r (y) 
(which are related by the definition (3.16)) as 

H { p\a 1 bi,a 2 b 2 , - ■ ■ ,a p b p ) = ((u r ) aibl (u r ) a2b2 ■ ■ ■ (u r ) apbp ^ (3.19a) 

/^(ai&i, a 2 b 2 , ■ ■ ■ , a p b p ) = {(v r )aib 1 (v r )a 2 b 2 • ' • (^a^y- (3.19b) 
We assume, for simplicity, that all odd moments vanish, i.e., 

/4t+i(ai&i, a 2 b 2 , • ■ • , a 2t+ ib 2t+1 ) = 0. (3.20) 
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FIG. 2: Schematic representation of the construction of the building block (BB) as a collection of 
"thin potential slices". 



We thus have 



( 



((u r )ab) = IH (ab) = 
u r )ab(us)cd) = ^ (ab, cd) 5, 



rs 



(3.21b) 



(3.21a) 



and similarly for the v r 's. It is useful to introduce the correlation coefficient between the 
matrix elements (u r ) a b and (u r )cd, (which coincides with the correlation coefficient between 
(v r ) ab and (v r )cd) as 



where $ (ab) = fi^ (ab,ab) denotes the variance of (v r )ab (recall that n " (ab) = 0). For 
even moments higher than the second we do not make, at this point, any special assumption; 
a particular scaling law will be assumed in Eq. (3.39) below. 

From the statistics of the (ii r )abS (and (v r ) a bs) we can find the statistics of the (e r )* 3 b , 
using the relations (3.14). For instance, we find that the first moment of (e r )^ 6 vanishes, i.e., 



C(ab, cd) 



/4 (ab, cd) /4 (ab, cd) 



(3.22) 



[/4 W) (ab)^ (cd)] 1/2 (ab)^ (cd)] V2 ' 




(3.23) 
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and that the second moments can be written as 

((eryi(e s ) h S) = l£\ab,cd) <U (3.24) 

where (i?r)£j was defined in Eqs. (3.14) and (3.15). The individual transfer matrices depend 
on the slice position x r and, as a consequence, they are not identically distributed. 



2. The transfer matrix for the Building Block. Its first and second moments 
The transfer matrix for the total sequence of m delta slices is given by 

M (m) = M m M m ^ l ---M 1 (3.25a) 
= (J 2i v + e m ) (I 2N + e m _i) • ■ ■ (J 2 jv + ei) (3.25b) 

= I 2 n + ^ ^ e r + ^ e n e r2 + ■ ■ ' 

•••+ e n---£r, + --- (3.25c) 

ri>...>r M 

= I 2N + e. (3.25d) 

The last line defines the matrix e [that was already introduced in Eq. (2.4)] by which the 
total transfer matrix M of the BB differs from the unit matrix I 2 n] it is given by 

e = ^2 e r + '^2 Sr i €r 2 1" e n'" £ rJ • (3.26a) 

r r\>T2 ri>...>r M 

m 

= J>W (3.26b) 

where the last line defines the contribution to e of order \i in the individual e r 's. Our aim is 
to find the statistical properties -in particular the moments- of the matrix e. In the future 
we shall use the notation (■ ■ ■ ) SL to indicate an average associated with the BB, i.e., 

(G(M))sl^(G(M^)), (3.27) 

just as in Eq. (2.8). For the average of M we trivially find, from Eqs. (3.25a), (3.25b) and 
the fact the various e r ' are statistically independent and average to zero [Eq. (3.23)], 

(M) SL = (M m ) ■ ■ ■ (M 1 ) = I 2N . (3.28) 
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Thus Eq. (3.25d) implies that the first moment of e vanishes, i.e., 



SL 



0. 



as could also have been obtained by averaging Eq. (3.26) directly: 

m 

<«>« = E 



'SL 



E <*> + E < 

) + •••+ ^2 ( e n"" e r P ) + 



ri>r 2 



ri>...>r M 



(3.29) 



(3.30a) 



(3.30b) 



For the second moments of e we have, from Eq. (3.26b) 



\ t ab t cd/ 



SL 



r 



E < l^} 

(1)1 a \jiv m \ 



(/*') 



id 



cd 



SL 



SL 



rare) 

+ < re [^'Q SL + (re + (re t^i 



<5L 



51 



(3.31a) 
(3.31b) 
(3.31c) 



+ 



The second line, Eq. (3.31a), is second order in the individual [e r ]^ b and hence in the 
potentials (v r ) a b, and the successive lines are higher order in these quantities. 

a. The second-order term in the second-moment expansion, Eq. (3.31a) The second- 
order term, Eq. (3.31a), in the second moment expansion can be written using Eqs. (3.26) 
and (3.24) as 



I cd 



SL 



r,s 
r,s 



(3.32a) 
(3.32b) 

(3.32c) 



From the definition of the correlation coefficient between pairs of matrix elements, Eq. 
(3.22), we can write the fraction in Eq. (3.32c) as 



{oh, cd) 
d 



C(ab, cd) 



$\ab) $\cd) 



d 



1/2 



C(ab, cd) 
V^ab(k) l cd {k) ' 



(3.33) 
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Here we have used the standard definition of the mean free path (mfp) £ a b associated with 
the incoherent sum of reflections from channel b to a from a sequence of v — 1/d scatterers 
per unit length, i.e., 

S^HlM*))-!'). (3-34) 
together with the fact that the average reflection coefficient for a delta slice is r-independent 
and approximately given, in the weak-scattering regime, Eq. (3.17), by [see Eqs. (A8), 
(3.13) and (3.14)] 

(\(riU\ 2 ) « ([(£iU] 2 }. (3.35) 
We can write the following equivalent expressions for the inverse mfp: 



(v) f , x _ ( ab ) _ /4 M) ( Q ft) _ /4 U) ( ab ) 

ujj.2 [ao) — - — — — - — _ — — - . 



(3.36) 



£ab(k) d Ak a k h d Ak a k b 

where the energy dependence of the mfp is exhibited explicitly. In the last member of the 
above equation we have defined the quantity pL^i^ab) as 

ft\ab) = . (3.37) 

a 

Since our delta slice is spatially symmetric in the x direction, we have the same result for 
the mfp for the transmission, out of the incident flux, from channel b to channel a. Within 
the present model there is thus no distinction between the so called transport and scattering 
mfp 's [26]. 

We now turn to the summation in Eq. (3.32c). We shall evaluate it in the dense-weak- 
scattering limit (DWSL) which we now define (see Eqs. (3.38) below). This limit was already 
referred to in Sees. I and II. Within the regime defined by the inequalities (3.1a) we have 
already considered a as the smallest length scale occurring in the problem and simplified 
the situation by literally taking the limit a — > [Eq. (3.8b)]. With regards to the next 
length scale in our regime, i.e., the distance d between successive scattering slices, we shall 
again be interested in a simplifying limit. For a fixed energy (and hence fixed X), fixed 5L 
and mfp's, it will be convenient to take the continuous limit 

d -> (3.38a) 

m —>■ oo, (3.38b) 

in such a way that 

md = 5L (3.38c) 
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remains fixed. In order to keep the mfp's £ a b of Eq. (3.36) finite, we carry to an extreme 
the weak-scattering condition (3.17) and, for a fixed energy, take the limit in which the 
individual scattering units become infinitely weak, i.e., 

fi { 2 u) (ab) -> 0, (3.38d) 

in such a way that, together with d — > 0, Eq. (3.38a), the ratio jl 2 u \ab) of Eq. (3.37), and 
hence the i a b{k) of Eq. (3.36), remain fixed. This limit has to be considered as the extreme 
idealization of the inequality (3.17) and of the inequality d <C {X,5L,£} of (3.1a) for fixed 
energy, 5L and mfp's. 

We have already assumed in Eq. (3.20) that all the odd moments of u and v vanish. We 
shall now assume that the even moments, Eq. (3.19) with p = 2t, scale with d as 

/4" Vi&i, • ■ • , a 2t b 2t ) = d l /4tVi & i; ■ ■ ) a 2 ^2i), (3.39) 

/4" (°i&i> • • • ; °2t&2*) being independent of d, with a similar expression for 
/4^(ai&i, • • • , a, 2t b 2t ). Eq. (3.37) is the particular case of this last equation for t — 1. 

In the DWSL, the appearing in Eq. (3.32c) tends to an integral, which we denote by 



^£ d (K6L) ee hm £ [(# r )Z(#r)S ■ d 

r 

-8L/2 



$Z(x)$ h cd (x)dx , (3.40) 

SL/2 

where fi^ix) is given by Eq. (3.15) with x r replaced by x. We find explicitly 



AS>. SL) = (-)■+'■+' » , (3.41) 



K 



2 

a quantity with dimensions of length, K l ^ cd being given by 

«$& = + + (-l) h K + (-l) 1 * 1 ** (3.42) 
From Eq. (3.42), and using the notation of Eq. (All), we readily find the symmetry relations 

j^ij,hl _ Tfhl,ij _ _T^ij,hT /o ao\ 

so that 

^O*, 6L) = A h ^ b (k, 5L) = A|2(fc, 5L) . (3.44) 
22 



We thus have, for the expression (3.32) in the DWSL: 

a result valid for arbitrary k and SL. 

For the application to Eq. (2.9) we shall need the expansion of the moments of e in 
powers of SL, with the BB translated to the interval (L,L + SL); this will be done in Sec. 
Ill C below. For the time being we perform that expansion, for simplicity, with the BB 
centered at the origin. We see from Eq. (3.41) that the leading term of A^^(A;, SL) in an 
expansion in powers of SL is linear in SL [as is obvious from the integral definition itself, 
Eq. (3.40)], i.e., 

SL) = (-) i+h+1 5L + 0{SLf. (3.46) 

As a result, Eq. (3.45) shows that the leading term in an expansion in powers of 5L of the 
second-order contribution to the second moments of s for the BB behaves, in the DWSL, as 

lim (\eU] ij h [ e «l"\ = (-y +h+1 C{ab ' C £U5L + Q(5L) 2 . (3.47) 

DWS \l labl ic d / SL ^UW)Uk) 

b. The fourth-order term in the second-moment expansion, Eq. (3.31c) A similar anal- 
ysis is performed in App. C, Eq. (C2), for the fourth-order contribution to the second mo- 
ments of e, Eq. (3.31c): it is shown that the leading term of such a quantity, in an expansion 
in powers of SL, behaves, in the DWSL, as (SL/£) 2 , where i denotes a typical mfp [see Eq. 
(C6)]. From this result and Eq. (3.47) we thus have 

Jj c hi\ _ i \i+h+i C{ab, cd) 

The analysis of the two above particular cases is generalized to arbitrary moments in 
App. D. For an even moment (p = 2t) in the DWSL, the lowest-order term in Eq. (Dl) 
(this term is of order 2t in the u r 's) has a leading term in an expansion in powers of SL 
which behaves as {SL/t} 1 . Higher-order terms in (Dl) are higher order in SL. Also, the 
dependence on the cumulants of the potential higher than the second drops out in the DWSL. 
The contribution to the second moments obtained above, Eq. (3.48), represents, for t = 1, 
a particular case of this general result. For an odd moment (p — 2t + 1), the corresponding 
term behaves as (SL/£) t+1 . 

In conclusion, this proves the behavior of the moments of e that was mentioned in Sec. 
II, right after Eq. (2.9). 
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= (-r +h+1 ^== 6L + °( 6L ) 2 ■ ( 3 - 48 ) 



C. The Diffusion Coefficients and the Diffusion Equation 



We now generalize the above analysis to the situation in which the BB lies in the interval 
(L,L + SL). The integral in Eq. (3.40) has to be performed in that interval [the notation 
(• • • )l,sl in Eq- (2.9) and in some of the following equations indicates this fact] and Eq. 
(3.45) becomes 

lim (\ e M]\ \e^] h \) = ° {ab ' ° d) A^(k, SL)e iK ^ L+ "\ (3.49) 

DWSV iabl ^ d / L ,5L ^/ £ ah (k)£ cd (k) ^ 

while the expansion in Eq. (3.48) [taking into account Eq. (C8)] is now 

ij,hl 



DWS v "~ ~"'L,8L 

+ 



(SL) 



2 



+0(5L) 3 , (3.50) 

where K a ^ d was defined in Eq. (3.42). In Eq. (3.50) we have defined the "diffusion 
coefficients" D^(k, L) 

D ab h cd( k > L ) = H i+h+1 ^ =e^'> , (3.51) 
ab ' cd 2y/£ ab (k)£ cd (k) 

which depend on the energy (through the energy dependence of the mfp's and through 
^abcd) an d a l so on the length L. Notice that the diffusion coefficients are, in general, 
complex numbers; this, however, should not worry the reader, because the evolution of real 
observables will always turn out to be real [see, for instance, Eq. (4.5a) below]. 

From the relations (3.43) we readily find for the diffusion coefficients the symmetry prop- 
erties 

D%%(k, L) = D%%(k, L) = [D^ b %(k, L)}*. (3.52) 
We introduce the expansion (3.50) and a similar one for higher moments of e on the r.h.s. 
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of Eq. (2.9), thus obtaining a power series in SL: 

(F{M)) L+5L = (F(M)) L + {^(k,L)SL 



ijhl,Xn 
abed, a/3 



+ 



a'/3',AV' 
8 2 F(M) \ 



(SLf 



p dM; 



+0(SL) 2 . (3.53) 

The curly bracket in this last equation corresponds to the BB second moment of Eq. (3.50); 
the contribution [which starts with (SL) 2 ] of the third and higher moments is just indicated 
in the last line. We also perform on the l.h.s. of Eq. (2.9) a Taylor expansion of (F(M)) l+sl 
in powers of SL around the "initial" value (F(M))l, i.e., 

{ F(M ))ml = (F(M)) L + SL + m> + . . . . (3.54, 

We then identify the coefficients of the various powers of 5L in Eqs. (3.53) and (3.54). In 
particular, the coefficients of SL give the diffusion equation, Eq. (2.10), derived in Sec. II, 
which we reproduce here: 

tjhlXfj, \ act cp I l 

abcdctf3 

Equating the coefficients of higher powers of SL in Eqs. (3.53) and (3.54) we obtain results 
which could be derived from the diffusion equation (3.55) by successive differentiations. We 
have verified this statement explicitly for the coefficients of {SL) 2 in the specific one-channel 
case treated in Sec. IV A below. 

The diffusion equation (3.55) governs the evolution with length of the expectation value 
of physical observables. The expectation values appearing in Eq. (3.55) must fulfill, for 
L = 0, the "initial condition" 

(F(M)) L=0 = F(I), (3.56) 

obtained by setting M = I in the expression F(M) for the observable, since for L = the 
scattering system is absent. More general initial conditions are discussed in Ref. [6]. 
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As was indicated earlier, the cumulants of the potential higher than the second are irrele- 
vant in the end; this signals the existence of a generalized central-limit theorem (CLT): once 
the mfp's are specified, the limiting equation (3.55) is universal, i.e., independent of other 
details of the microscopic statistics. 

Since the structure of the present diffusion equation is essentially the same as the structure 
of the one derived in MT (Ref. [20], Eq. (3.18)), it is worthwhile, for the sake of comparison, 
to summarize, at this point, the MT model. In MT the statistical assumptions are made 
at the level of the individual scattering units, just as in the present paper (the same units 
that were also contemplated in Ref. [12]); however, the assumptions are not made for the 
potentials, but rather for the corresponding transfer matrices. In MT, the transfer matrix 
for each scattering unit is close to the unit matrix and is written as M r = I + e r , just as 
in our Eq. (3.13) above; it is further expressed in terms of independent parameters (in the 
Pereyra representation [27]), for which various statistical assumptions are made: 

i) The first moment and some of the second moments of the independent parameters are 
chosen so that the resulting (e r ) = (see Eqs. (3.15) and (3.16) of MT; with this feature, 
there is no drift term in the resulting Fokker-Planck equation), while the remaining second 
moments of the independent parameters are kept arbitrary. 

ii) The individual scattering units are statistically independent and identically dis- 
tributed. 

iii) The energy does not appear explicitly, but only as the energy at which the resulting 
mfp's have to be evaluated. 

iv) In order to obtain explicit expressions for the diffusion coefficients, in the analysis that 
follows from Eq. (3.18) of Ref. [20] a more explicit model was postulated for the second 
moments mentioned in i) above. 

In the present paper, assumption i) is a consequence of the vanishing of the first moment 
of the individual potentials, Eq. (3.21a), thus giving Eq. (3.23). Assumption ii) has to 
be contrasted with Eq. (3.24) above, which shows that, here, the transfer matrices for the 
individual scattering units are not identically distributed. As it has already been stressed, 
in contrast to assumption iii) the energy appears now explicitly. Finally, the additional 
assumptions mentioned in iv) are, to some extent, arbitrary; they are compared below with 
those arising from the short-wavelength approximation of the present model. 
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D. The Short-Wavelength Approximation. The regime (3.1b). 



In the DWSL the above expressions are exact for all energies. We now turn to a different 
regime, to be called the short-wavelength approximation (SWLA), defined by the inequalities 
(3.1b). The regime to be studied is analogous to the geometrical optics limit studied in optics 
[28]. Essentially, we shall assume that we can fit many wavelengths inside a BB, i.e., 

A < 5L , or kSL > 1 , (3.57) 

so that in this regime only lengths much larger than the wavelength actually enter the de- 
scription. 

To this end we go back to Eq. (2.9) which, after setting {^ a k c ) LSL = because of (3.29), 
we rewrite here for convenience: 

(F(M)) L+SLik = (F(M)) L>k + i E v4 e%) w E M % Zl^X ) 

' ijhl X/j, \ 0m aa ° m c(3 / L>k 

abed ctf3 

+ •••. (3.58) 

We have indicated explicitly the k dependence of the various expectation values. We first 
analyze below the BB factors appearing on the r.h.s. of the above equation, and then the 
remaining expectation values. 

1) The BB factor (e l ^ h E%d) LSL . k can be written, from Eq. (3.31), as 

<4 «S> w = mi TO W + W) W + (3-59) 

The first term on the r.h.s. of this last equation is given by Eq. (3.49), and its contribution 
to (3.58) is given by 



ijhl X/i 
abed a/3 



1 C(ab,cd) Aij,M(u xt\Jk(l+^)\^ /( \ijhl\p v 

2 ijhi vM*0M*0 v 

afecd a/3 
(tf=0) 



1 C(ab,cd) ij:M . xt\Jk(l+^)\^ /( yjhix»s 

abed ct/3 



ijhl V -aow-cavv A/x 

(3.60) 
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In this equation, K is an abbreviation for K a ^ d which was defined in Eq. (3.42), and 
A^'^fc, SL) was given in Eq. (3.41). We have also used the notation 

• ywA " \ = / M jX M l " 92 F ^ \ (3 si) 

\ Ulvl aa ulvl cl3 I L k 

as an abbreviation for the last factor appearing on the r.h.s. of Eq. (3.58). In the one but 
last line in Eq. (3.60) the sum is over the combinations of indices that make K a ^ h cd = 0, 
while in the last line it is over those combinations that make K a ^ d 7^ 0. 

The second term on the r.h.s. of Eq. (3.59) is given in Eq. (C7) and, using a similar 
convention as in the last equation, its contribution to (3.58) can be written as 

ijhl A/i 
abed a/3 



1 C(aa\ cp) Cja'b, fid) 

2 Jfa ^£aa<(k)e c p,(k) VWfc)W*0 



abeda' P' 
(K 1= K2=0) 



AiX',X'j,hfi',n'l / 7 -r>/ rn\ i(K 1 +K 2 )(L+^) \^ I ( \ijhl\p \ 



A/ 1 

a/3 



+ lj ^ + + ^ I C(aa',c(3) C(a%(3'd) 



2 I ijW Wy! Wyi J V^AWcAkl vWP) 

abeda' P' abeda' '/?' abeda' f3' 

(^17^0,^27^0) (^1=0,^27^0) (ATi^O,^^) 



X 



7e(5L))e^+W-^> £<(• • OjfflW • (3-62) 

A/; 

We recall that K\ and are defined in Eq. (C4). 

Higher-order contributions occurring on the r.h.s of Eq. (3.59) can be obtained from the 
analysis of App. D. 

We now analyze the consequences of the inequality (3.57) for the above expressions (3.60) 
and (3.62), which so far are exact. It will be convenient to take the wavenumber k as 

k = (5^, (3.63) 
i.e., halfway between the threshold for the last open channel and that for the first closed one, 



so that the longitudinal momenta are given by k a = ky/1 - [a/{N + 1/2)] 2 . From Eq. (3.42) 
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we see that when K l J b ^ d ^ 0, K l J^ is proportional to k (the coefficients only depending on 
channel indices), so that K l J b ,h cd 5L 1. As a result: 

i) in Eq. (3.60) the sum with K = gives the largest contribution (proportional to 5L, as 
we now analyze in detail), while the sum with K ^ 0, which contains K in the denominator 
of A^ 6 ' crf (A;, 5L), will be neglected. 

Let us be more specific about the combination of indices ab, cd and ij, hi that give rise 
to K = in Eq. (3.60). Take, for instance, i = j — h — I = 1. Since k a , kb, k c , kd are 
incommensurate, K^'^ = kb — k a + kd — k c (see Eq. (3.42)) can only vanish if a = b and 
c = d, or a = d and b = c. On the other hand, K^'^ = —(k a + kb + k c + kd) never 
vanishes. We thus have, for A^ fe '^(/c, 5L), defined for arbitrary k and 5L in Eq. (3.41), the 
approximate result: 

8L) « (-iy +h+1 5 K0 ■ 5L, (3.64a) 

(here, 5kq is Kronecker's delta which takes on the value 1 when K — and vanishes other- 
wise) or, more explicitly: 

A 11,11 / / XT \ _ Sabred + 5 a d5 be . . 

A «6,«i(M£) ~ 7—? 5L ( 3 - 64b ) 

a 11,22 / ; r r \ dabbed + "ac^ta xr (naA \ 

A abcd( k , 5L ) ~ 7— T 6L ( 3 - 64c ) 

J- + Oad 

AgS(ML) « ^m**** 6 * 5L (3.64d) 

«<) « 4£;2(ML) « Ag2(*, *L) « . (3.64e) 

We can thus write / [e {1) ]^ b in the DWSL, followed by the SWLA, as 



lim / \e<»]\ \e^] k [) » g^fg (-l) <+h+1 <fco • *L, (3.65a) 
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One finds explicitly in the various cases (C(a, c) being an abbreviation for C(aa,cc)): 
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I cd 
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1 + 5 ad 



C(a, c 
C(a, c 



, n hO rf j H " — OndMhr 



5L 



1 + 5, 



ab ^ab 



(3.65b) 
(3.65c) 
(3.65d) 
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lim (le^Ylle^} 1 ') « 0. 
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(3.65e) 



where we have defined the diffusion coefficients in the SWLA as 



jy,j,hl 
ab,cd 



1 y+/ t +i C(ab,cd) g 



2^£ ab (k)£ cd (k) 



AT) 



Other combinations can be found from TRI, Eqs. (A10). The result is that in the DWSL, 
followed by the SWLA, the second-order contribution to a second moment of e for the BB is 
either negligible or behaves as 5L/£, £ denoting a typical mfp. 
One can write Eqs. (3.65) as 

(3.66) 



(3.67) 

(3.68a) 
(3.68b) 
(3.68c) 
(3.68d) 



which, from Eq. (3.65), take the explicit form 



^ab\cd (k) 

Dltfd(k) 



1 + 8 n 



1 



OabOcd + -^—OadObc 



2y/£ aa (k)£ cc {k) 

C(a, c) 



2(\ 



ab 



1 



=5 n h5 rr i -\- — - — 5 nr 5hd 

2^£ aa (k)Uk) 2£ ab 



l + 5ad 
0~acO~bd + 0~adO~bc 1 

l + 5ab 2£ ab (k) 



,11,21 



12,12, 



These diffusion coefficients depend on the energy through the mfp's only. 

ii) Eq. (C3) shows that in the DWSL, followed by the SWLA, the fourth-order contribution 
(C2) to a second moment of e for the BB is either negligible or behaves as (SL/£) 2 , £ denoting 



30 



a typical mfp. Thus in Eq. (3.62) we keep only the sum for K\ = K2 = and neglect the 
other summations, the result being thus proportional to (5L/£) 2 . 

We finally obtain, for the BB second moments of Eq. (3.59) in the SWLA: 



a'/3',X'iJ,' 

+0{5Lf . (3.69) 



2) Similar arguments applied to the analysis of App. D lead to the result that a (2t)-th 
moment of e for the BB can either be neglected because it contains fc a 's in the denominator, 
or it gives a contribution to Eq. (3.58) which is porportional to (SL/£Y, whereas a (2i + l)-th 
moment contributes as (SL/£) t+1 . 

3) We need some knowledge about the behavior of the averages (• • - )h,k appearing in Eq. 
(3.58) in the SWLA. We shall assume that, for large enough k, we can approximate 

(■■■>L,^(--->i 0) , (3.70) 

where the r.h.s. represents a function smooth to all scales of L and whose energy dependence 
only appears through the mfp's £ a b(k). This ansatz, which seems merely reasonable at this 
point, is verified in a particular case in Sec. IV A below. In the analysis that follows we 
shall assume that the energy is kept fixed, so that the mfp 's will be taken as fixed parameters 
and will be written as £ a t,. Likewise, we shall write D l ^ d for the diffusion coefficients. 

We now make use of the results li), lii) and 2) above, as well as the assumption (3.70), 
to write Eq. (3.58) in the SWLA as 

(°) ~ /cVA/n\(°) 



+ 



E 

ijhl.Xfi 

abcd,af3 



Mi 



XM[ &F(M) \ (0) 



lba ^dM£dM%/ L 
+0(5 L) 2 . (3.71) 

The square bracket in this last equation corresponds to the BB second moment appearing 
in Eq. (3.58); the contribution of the third and higher moments is just indicated in the last 
line, in accordance with 2) above. 
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We now assume that the quantity {F{M))f ) +SL appearing on the l.h.s. of Eq. (3.71) can 
be expanded in a Taylor series around the value L, and that 5L is smaller than the radius 
of convergence R of the expansion, i.e., A <C 5L < R, so that: 

(F(M))fl SL = <F(M)><°> + 8 -Wll SL + ^ {F ^ ))f ' (SL)' + .... (3.72) 
Comparing the coefficients of 5L in Eqs. (3.71) and (3.72) we finally find: 

abcda/3 

In the SWLA we have thus ended up with an evolution equation for the "smooth" quantities 
defined in Eq. (3.70). 

We need to fix the initial conditions appropriate to Eq. (3.73). If we require Eq. (3.56) 
for the exact expectation values, i.e., (F(M)) L=0>k = F(T), and F(I) is fc-independent, then 
Eq. (3.70) implies 

(F(M))i°! = F(I). (3.74) 

More detailed assumptions than (3.70) on the structure of the expectation value (■ ■ - )L,k 
appearing in Eq. (3.58) in the SWLA are presented in App. F for the one-channel case, 
N — 1. There, a rederivation of Eq. (3.73) using such assumptions is also discussed. 

The derivations given above of both diffusion equations, Eq. (3.55), valid for arbitrary 
energies, and Eq. (3.73), valid in the SWLA, use, as a starting point, Eq. (3.58), which 
describes the result of adding a BB to an already existing waveguide of length L. This is 
also the starting point of the derivation given in App. Fl. We believe that it would be 
very instructive to rederive the diffusion equation in the SWLA, Eq. (3.73), starting directly 
from the more general one, Eq. (3.55), since such a derivation would shed more light on the 
nature of the various approximations involved. However, we have succeeded in fulfilling this 
goal only in the one-channel case, N = 1; the derivation is presented in App. F2. 

IV. APPLICATIONS OF THE DIFFUSION EQUATION 
A. Analytic examples 

In this section we study a simple example in which the diffusion equation (3.55) can be 
solved exactly. We restrict the analysis to a one-channel geometry (N = 1) and consider, as 
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examples of the observable F(M), the quantities 

1 1 



M n M 22 = aa* = — = - (4.1a) 
tt* T 



M n M 12 = a /3=-(^j , (4.1b) 

where we have used Eq. (A8) to establish the connection with reflection and transmission 
amplitudes. We shall give only the main results of the calculation, some of the details being 
presented in App. E. 

For the one-channel case, the diffusion equation (3.55) can be written as 

= V n L) /*» M lti »JM ) , (4.2) 

ijhlX/j, x 1 L 

where the diffusion coefficient D^ ,hl (k, L) is given explicitly in Eq. (El). For simplicity, we 
have suppressed all channel indeces, which would take the value 1. We emphasize that in 
the DWSL this equation is exact, in the sense that it is valid for all energies. 

The mfp is energy dependent. However, in the present calculation we keep the energy 
fixed and so the mfp is taken as a fixed parameter and will be written as £. One can write 
all the evolution equations in terms of the ratio of the length L to the mfp £ 

s = L/£ , (4.3) 

and essentially the ratio of the mfp to the wavelength A 

x = 2k£. (4.4) 

Using the diffusion coefficients of Eq. (El) one finds the pair of coupled equations 

d(aa*) s 



ds 



(a(3) s e lX0S + (2 {aa*) s - 1) + (a*/3*) s e~ tx ° s , (4.5a) 



= - (aP) s ~ (2 (aa*) 8 - 1) e~ lX0S - (a*P*) s e~ 2iX ° s , (4.5b) 

which have to be solved with the initial conditions at s = 0: 

(aa*) s=0 = 1 (4.6a) 

(aP) s=0 = 0. (4.6b) 



33 



The second derivatives of the observable F(M) appearing on the r.h.s of the diffusion 
equation (4.2) produce, in general, quantities which are different from the observable F(M) 
itself, whose average we wish to study. One then needs to compute the evolution of these 
other quantities and this, in turn, generates still new ones. In the example considered here, 
Eq. (4.5) shows that the evolution of (««*) involves (aa*) and (a/3), and similarly for the 
evolution of (a/3): we thus find a pair of coupled equations which "close", in the sense that 
the quantities occurring on the r.h.s. are the same as on the l.h.s. 

The evolution equations (4.5) for the real quantity (aa*) s and the complex quantity 
{a/3) s can be written as the triplet of coupled equations (E2), which can be solved using the 
method of Laplace transforms, using the initial conditions (4.6), with the result 



(«« >, = 2 + 2 



Pi + 2pi + x 2 Q cPiS p\ + 2p 2 + x 2 cP2S p\ + 2p 3 + x\ cP3S 



. (Pi - P2) (Pi - Ps) (P2 - Pi) (pa - Ps) (P3 - Pi) (Pa - P2) 

(4.7a) 

Pi + lX e (pi-ix )s + Pi + ZX c (P2-ix )s 



(Pi - P2) (Pi - Pa) (P2 - Pi) (P2 - Ps) 



+ ^0 c ( P3 -ixo) 



(4.7b) 

(P3 - Pi) (P3 - P2) 

In this equation, pi, P2 and P3 are the roots of the third degree polynomial P (p) = p 3 + 
XqP — 2xq, with pi e R, P2,P3 e C and p% = p* 2 . 

The solutions (4.7) are exact, being valid for arbitrary length L, mfp I and wavenumber 
k. Moreover, as shown below, the solutions of the diffusion equation are in full quan- 
titative agreement with the statistical averages obtained from numerical solutions of the 
one-dimensional wave equation. 

A one-dimensional version of the delta-slice model discussed in Sec. Ill B 1 is sketched 
in the inset of Fig. 3. Notice that in a ID problem there are no evanescent modes. The 
system of length L is constructed from "delta potentials", U r (x) = u r 5(x — x r ) (recall that 
U r (x) and u r have dimensions of k 2 and k, respectively), located at the positions x r = rd 
(r = 0, 1, 2 . . . ); u r is assumed to be uniformly distributed over the interval [— u max , +u r 
The mean free path, obtained from Eq. (3.36), is simply given by: 

d 1 / v \ 2 



'"max] 



(4. 

I 3 V 2k ) 1 
The results of the numerical calculations for (aa*) and (af3) versus L are shown in Figs. 3 
and 4, respectively. Averages were obtained from 10 7 different microscopic realizations. Nu- 
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100 



FIG. 3: (M 11 M 22 ) = (aa*) versus kL. Numerical results (bold line) from the one-dimensional 
model sketched in the inset are indistinguishable from the analytical results, Eq. (4.7a). The 
results correspond to xq = 2k£ = 200 and d/i = 10 -3 . 



merical results (bold line) are indistinguishable from the analytical solution of the diffusion 
equation (Eqs. (4.7a) and (4.7b)). 

It will be interesting to see what these results reduce to in the SWLA discussed in Sec. 
HID above. In preparation for this, we first consider a fixed value of s = L/i and take 
xo = 2k£ ^> 1. From Eq. (E8) one can expand the functions (aa*) s , (a(3) s in powers of 
1/xq; in terms of the original variables k, I and L they take the form 



(aa*) , 



(2kiy 



L\ l l e 2ikL + e~ 2ikL 
1 + 2— ] e 2 ~ + e~~ 



(4.9a) 



e «e 



-2ikL\ 



) (2key 



5-3^ 



_L 

■e * 



e i e 



-2ikL 



-O 



k(: 



(4.9b) 



The solutions (4.9) satisfy the differential equations (4.5) together with the initial conditions 
(4.6) to every order in the expansion in powers of l/2k£. 

Notice that the ansatz made in Eq. (3.70) is verified explicitly in this example, with the 
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FIG. 4: Real (top) and imaginary (bottom) parts of (M n M 12 ) = (a{3) as a function of kL. 
Numerical results (bold line) are indistinguishable from the analytical results. The parameters are 
the same as in Fig. 3. 

result 

<««*>i 0) = ^(1 + e 2 ^), (4.10a) 
( a p)f> = 0, (4.10b) 

which represents, in this particular case, the SWLA discussed in Sec. HID. The result 
(4.10a) agrees with what had been obtained earlier as a solution of the diffusion equation of 
Ref. [15] for N — 1, also known as Melnikov's equation. Notice that 

/D\(0) 1 

</Wf = ("Of - 1 = (|) = - 1) (4.11) 
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represents the well known exponential increase of Landauer's resistance [29] . 

If, in Eq. (4.9), we further expand the exponentials e 2L ^ e , e _i//£ , in powers of l/£, we end 
up with an expansion of (aa*) s , (a/3) s in powers of l/£. The result found in Eq. (3.49), 
setting L = and interpreting 5L as L, is precisely the term proportional to l/£ in such an 
expansion; we have verified the consistency of the two results up to 0(1/1). 

B. Random walk in the transfer matrix space: Numerical simulations. 

As we have shown, the diffusion equation, Eq. (3.55), determines the statistical properties 
of transport for any physical observable and it only depends on the mean free paths £ a b. Once 
the various £ a b are specified, the statistical distributions are universal, i.e., independent of 
other details of the microscopic statistics. However, in order to know the exact shape of 
the distribution of a given observable we have to solve the diffusion equation. This is a 
challenging problem even in the isotropic case [15] (where all the mfp's are equivalent, 
£ab = £)■ Here, instead of a direct solution of the multidimensional diffusion equation we 
have followed an alternative way that can be seen as a generalization of a random walk 
in the transfer matrix space (Ref. [30]). The method, based on our previous theoretical 
description, can be summarized as follows: 

1) We first obtain a set of mean free paths from a given microscopic potential model for 
the building block or, eventually, from specific experiments on very thin slabs. 

2) We generate an ensemble of transfer matrices having their first and second moments 
equal to those corresponding to a BB of a certain length 5L. 

3) The transfer matrix for a system of length L = PSL is obtained by combining P 
building block matrices randomly chosen from the ensemble. This procedure can be repeated 
again and again in order to obtain the statistical distribution of any physical quantity. As 
predicted by the CLT associated with the composition of BB's explained in App. G, higher 
order moments of the BB matrix elements play no role in the final statistics. 

The statistical distributions of different physical quantities will be shown to be in full 
agreement with the results of exact microscopic numerical calculations for a model system. 
This shows that validity of the diffusion equation given in Eq. (3.55) goes beyond the various 
formal limits discussed in Sec. III. 
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1. Microscopic potential model and mean free paths 



Let us consider the potential model sketched in Fig. 5. In this model, a 2D waveguide 

8jc 



by 



— 



5 W 



W 



FIG. 5: Schematic representation of the microscopic model based on random potentials. Each 
square of the plot, or "cell", represents a region of constant random potential. 

with perfectly reflecting walls has a region of length L which is divided into small "cells" of 
dimensions 5x x Sy. The working wavelength is chosen to be such that Sx, 8y <C A. In the 
language of Sec. Ill A, the potential in the r-th slice, Eq. (3.9), is replaced here, for finite 
5x, by 

U r (x,y) = Ur(y) 9&e{X 5 ~ Xr \ (4-12) 

where 6$ x (x — x r ) takes the value 1 inside the interval (x r — Sx/2, x r + Sx/2) and outside. 
Should 5x tend to zero, the expression in (4.12) would tend to that of Eq. (3.9a). Inside 
the r-th slice, the potential is taken to be constant within each cell, i.e., 

u r(y) = 5^ U » *v(V ~ f»0, ( 4 - 13 ) 
s 

so that 

U r (x, y) = } j U s 9 Sx (x - x r )9 Sy (y - y r ), ( 4 -14) 

s 

with U s = u s / 8x. The constant values U s of the potential inside each cell located in the region 
W — 8W < y < W is sampled from a uniform distribution within the interval [—Uq,Uq]. 
Outside the region defined by W — 5W < y < W, the potential is taken to be zero. 

In order to get the mfp's corresponding to our model system, we follow the same steps 
leading to Eq. (3.36) in Sec. Ill above. In the limit 5x, 5y 5W, and neglecting the 
coupling to evanescent modes, i.e., using the "bare" potential u instead of the "effective" 
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one u (see text following Eq. (3.16) and App. B), we obtain 

1 ({v ab } 2 ) U 2 SxSy 



xi(y)xi(y)dy, (4.15) 

'-ab OX 6 4K a Kb JwSW 

where Xa{y) are the transverse eigenfunctions of the clean waveguide (Eq. (3.3)). The mfp's 
for bulk disordered systems, i.e., when the disordered potential covers the whole section of 
the waveguide (5W = W) are simply given by 



1 



U% Sx5y 2 + 5 ab 



bulk 3 4* a A* 2W 



(4.16) 



In order to analyze a surface disordered waveguide, we shall also consider the limit 5W <C W, 



1 



Ul 5x5y ^4vr 4 



surface 3 4k a k b V W 2 



2. Random transfer matrices for a building block 

In order to generate an ensemble of random transfer matrices whose first and second 
moments are given, it is useful to describe the transfer matrix elements of the BB as a 
function of the 2N 2 + N independent parameters of the Pereyra representation [see App. A, 
Eq. (A14)]. The matrix e of Eq. (3.25d) can be expressed (in that representation) as 

e 11 = e*V 1 + VV* ~ 1, (4.18a) 
e 12 = e ih r], (4.18b) 

where h is an arbitrary N x N Hermitian matrix (thus contributing N 2 parameters) and r) 
is an arbitrary N x N complex symmetric matrix (thus contributing N 2 + N parameters). 

Applying successive approximations to Eqs. (4.18) it is possible to invert them to express 
the matrices h and rj as functions of the blocks e % \ i.e., 

h = -is 11 + 1 (e 12 e 21 + e u e u ) + O (e 3 ) , (4.19a) 
n = e 12 - e ll e 12 + O (e 3 ) . (4.19b) 



The aim is to derive the statistical properties of the matrices h and rj in terms of those of 
the blocks e 1 ^ which we derived in the previous section; we shall do this in the SWLA (see 
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Sec. HID). We can use Eqs. (4.19), (3.68) and (3.69) to obtain (in powers of SL) the first 
and second moments of the matrix elements t] ab , h ab . For the first moments we obtain 



(hab) SL 
(Vab) SL 



and for the second moments 

(habh c d)s L = - (taWcd) SL + 

+0 (SL) 2 , 

(KbKd)dL = ~ ( £ ab £ fd) 6L + 

+0 (SL) 2 , 



)sL + 



(habVcd) SL = 


i ( £ al £ cd 


(habV*cd)sL = 


~i ( £ ab £ cd 


(VabVcd} SL = 


\ c ab° cd / SL 


(VabV*cd}6L = 


(f 12 f 21 \ 

X^ab^cd/SL 



O (SLf 

(SLf 

5L 

1 + Sac 

SL 



SabScd- 



C (aa, cc) S ad S bc 

's/^aa^cc ^ab 



s, 



ad 



C (aa, cc) S ac S bd 

OabOcd , „ „ + 



■ah 



■ = 0(SL) 2 , 

■ = 0(SL) 2 , 

O(SL) 2 , 

SacSbd + S ad S bc SL 



l + S 



ab ""ab 



O(SL) 2 . 



(4.20a) 
(4.20b) 



(4.21a) 

(4.21b) 

(4.21c) 
(4.21d) 
(4.21e) 
(4.21f) 



To generate the ensemble of random transfer matrices for r/ and h in the SWLA, we need 
to know the statistical properties of the real and imaginary parts of the matrix elements r] ao 
and hob, to be denoted as 

Vab = Rer/ a fe = - (r] ab + rf ab ) , rf ah = lmr] ab = — (r] ab - rf ab ) , (4.22a) 
h% = Reh ab = - (h ab + h* ab ) , hi b = lmh ab = - (h ab - h* ab ) . (4.22b) 
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Using Eqs. (4.20) and (4.21) we find 

((-fi)\ = <to) I ) tt = £ + < ,£ > > ' Va A (4.23a, 

(K) 2 ) eL = (KT) SL = ^ +0{SL)2 - a * b ' (423b) 

((^aaf) SL = t— + 0(<5L) 2 , (4.23c) 

(h aa h bb ) SL = C( ^!^ 6L + o (5Lf , (4.23d) 

V ^aa^bb 

(hX) gL = {h'jQn = {h R ab h{ d ) 5L = O (5L) 2 , a^b^c^d, (4.23e) 

W&L = (hMd) 8L = { h a b hL) 5L = O (5L) 2 , a^b^d, (4.23f) 

(hS,r&) SL = (KMsl = (habVDsL = (KbV^)s L = O (5L) 2 , (4.23g) 

{r&vDu, = { h abK d ) 5L = O (5L) 2 (4.23h) 

(r&r&) SL = (vLvDsl = O (SL) 2 , a + c, b + d, (4.23i) 

(VaMsL = W ad ) SL = (VaM SL = O (6L) 2 , a + b + d. (4.23j) 

We recall that the diagonal elements h aa are real since h is a Hermitian matrix. 

From now on, to generate the ensemble we shall consider a potential which is delta 
correlated in the transverse direction; in that case we have 

c ( QQ > hh ) = J_ ( 4 24 ) 

V ^aJbb tab 

which allows rewriting Eqs. (4.23c)-(4.23d) as one equation: 

(h aa h bb ) 5L = 5 -^ + 0(5L) 2 . (4.25) 

tail 

Therefore, in the SWLA, real and imaginary parts of the matrix elements of rj and off- 
diagonal matrix elements of h are, to order 8L, uncorrelated, with zero mean, Eq. (4.20), 
and with variance 8L/2£ ab , Eqs. (4.23a) and (4.23b). For these elements we have used two 
different distributions giving the same variances: 

1 



P^ x) 



Ay/3, 



a 



e(x + v^o) -e(x- V3a^) , (4.26a) 



P 2 (x) = U (x - a) + U (x + a), (4.26b) 

9 (x) being the usual step function and a 2 = var (x). As we can see from the CLT of App. 
G, the final results only depend on the coefficients proportional to SL, while the rest of the 
details of the distributions do not play any role. 
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In contrast, the diagonal elements of the h matrices are correlated, Eq. (4.25). In order 
to generate numerically a set of uncorrelated variables from the diagonal elements of the h 
matrices we have performed an orthogonal transformation on the diagonal terms h aa , 

b 

in such a way that the covariance matrix C a b = (h aa hbb) = o~L/£ a t> is diagonalized, to obtain 

(KaKb) = <W£ (4.28) 

Hence we can numerically generate a set of N uncorrelated variables h' aa with zero mean 
and a variance given by the eigenvalues of the C a b = 8L/£ a b matrix and, after that, obtain, 
by the change of coordinates (4.27), the h aa variables which are properly correlated. 

3. Random walk in the transfer-matrix space: Statistical conductance distributions 

Once we have numerically generated an ensemble of transfer matrices with its first and 
second moments correct up to order 5L, we can obtain a transfer matrix corresponding to 
a system of length L = P6L by multiplying P transfer matrices of the ensemble of BB's 
taken at random. Numerically this procedure is unstable because the pseudounitary group, 
to which the transfer matrices belong, is non-compact [10]. This property leads to numerical 
instabilities as the norm of the transfer matrix elements can grow without limit. Instead of 
using the product of transfer matrices, we obtain the scattering matrix associated with each 
transfer matrix (Eqs. (A8)), and then combine different scattering matrices to obtain the 
scattering matrix for the system of length L (Eqs. (A16)). 

For a given set of mean free paths £ a t, we choose the length 5L of the BB in such a 
way that SL/£ a b <C 1 for all channels. With this, we generate random transfer matrices as 
explained above and, for each one, we obtain the corresponding scattering matrix. Applying 
Eq. (A16) P times we obtain the scattering matrix corresponding to a system of length 
L = P5L. This procedure can be repeated as many times as needed to obtain the desired 
statistical properties. 

A detailed numerical analysis of the statistical properties is beyond the scope of the 
present work and will be discussed elsewhere. Here we just focus on the statistical distribu- 
tion of the conductance and the intriguing discrepancies between surface and bulk disordered 
systems [18, 19]. 
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FIG. 6: Bulk disordered waveguides: Average transmittances (Ta) (channel in = channel out), 
as a function of L/t. The inset shows the equivalent results for (Tij) with i ^ j. The results based 
on the numerical solution of the Schrodinger equation (microscopic calculation) and the random 
walk simulation of the diffusion equation in the SWLA overlap. 

a. Bulk Disorder The behavior of the average transmittances (Ta) (channel in = chan- 
nel out), for bulk disordered wires, is plotted in Fig. 6 as a function of L/£, i being the 
averaged transport mean free path, 



The inset shows the equivalent results for (Tij) with % ^ j. The random- walk simulation 
was performed in the SWLA. We have also solved numerically the Schrodinger equation 
for the same model system (sketched in Fig. 5). We followed an implementation of the 
so-called generalized scattering matrix (GSM) method (see for example Ref. [24]). The first 
step consists in the calculation of the set of transverse eigenfunctions and the scattering 
matrix for each slice of length 8x. The combination of two consecutive slices is done by 
mode matching at the interface. After that we combine scattering matrices to obtain the 
scattering matrix of the whole system. It is important to mention that this calculation 
is performed using both propagating and evanescent modes and hence, this method can be 
considered as exact. The statistical properties of any transport parameter obtained from 10 5 
different realizations were found to converge for 3 evanescent modes. The calculations have 
been done starting from the set of mean free paths £ a b given by Eq. (4.16) for kW = 5.5ir 
(corresponding to 5 propagating modes), UqW 2 = 100 and 5x/W = 5y/W = 1/50. The 




(4.29) 
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exact numerical results for the average transmission coefficients are indistinguishable from 
the random-walk simulations. 

The random- walk results for the average reflection coefficients (Rij) for bulk disorder 
(shown in Fig. 7) are also in good agreement with our numerical results as well as with 
previous numerical work [31] (using a two-dimensional Tight-Binding model with Ander- 
son disorder). The set of reflection coefficients corresponding to backscattering ((Ra)) are 
consistent with an enhanced backscattering factor (Ra) / (Rij) ~ 2, as expected from the 
DMPK equation. 



A 




FIG. 7: Bulk disordered waveguides: Average reflectance (Rij) as a function of L/l. The 
upper curves correspond to i = j and the lower curves to i ^ j. The results based on the numerical 
solution of the Schrodinger equation (microscopic calculation) and the random walk simulation of 
the diffusion equation in the SWLA overlap. 

The distribution of the dimensionless conductance, P(g) [with g = trace(ti')], for bulk 
disordered wires is plotted in Fig. 8 for different conductance averages, (g). The inset 
shows the average conductance as a function of L/i. The exact numerical results for the 
conductance distribution (histogram lines in Fig. 8) are indistinguishable from the random 
walk simulations (open circles). For comparison we also plot (continuous line) the exact 
result of the diffusion equation of Ref. [15] (DMPK equation) obtained from a Montecarlo 
simulation [18]. Despite the slight channel anisotropy of transport, the results are compatible 
with those of the DMPK equation. 
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FIG. 8: Bulk disordered waveguides: Distribution of the dimensionless conductance, P(g), 
for different conductance averages, (g). The three different curves based on different approaches 
overlap. Circles correspond to the random walk simulation of the diffusion equation in the SWLA. 
The continuous line represents the results of the Montecarlo simulation of Ref. [18]. The histogram 
lines are the results based on the numerical solution of the Schrodinger equation (microscopic 
calculation). The inset shows the average conductance as a function of Ljl. 

b. Surface Disorder In the case of surface disorder, the mean free paths are very dif- 
ferent from those obtained for a uniform (bulk) distribution of scatterers. In particular, the 
dependence of £ a b on a 2 b 2 [see Eq. (4.17)] reflects the strong channel anisotropy of transport 
in surface disordered waveguides [11, 32-36]. This could be the origin of the differences be- 
tween bulk and surface distributions. Previous numerical calculations for surface disordered 
waveguides, showed that close to the onset of localization, the conductance distributions 
presented an unexpected sharp cusp-like shape [19]. The distribution of the dimensionless 
conductance for surface disordered wires obtained from the random walk simulation in the 
SWLA is plotted in Fig. 9 (open circles) for different conductance averages. The exact 
solution of the Schrodinger equation (microscopic calculation; histograms) is again in full 
agreement with the diffusion equation and with previous numerical work [11, 19]. The cal- 
culations have been done starting from the set of mean free paths £ a b given by Eq. (4.17) 
for SW = 0.1W, U = 100/IF 2 , kW = 5.5vr (N = 5), 6x = 106y = IF/50. 

It is worth noticing that when the disordered region is confined close to the surface, the 
mean free paths can be extremely large (for example, for the present calculation, I pa 1.50 x 
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FIG. 9: Surface disordered waveguides: Distribution of the dimensionless conductance, P(g), 
for different conductance averages, (g). Circles correspond to the random walk simulation of 
the diffusion equation in the SWLA. The histogram lines are the results based on the numerical 
solution of the Schrodinger equation (microscopic calculation). The equivalent results for bulk 
disorder (continuous line, DMPK) are also shown for comparison. 

lO 4 !^). The exact numerical solution of the wave equation is then extremely expensive in 
terms of computation time compared to the random walk simulations based on the statistical 
properties of the BB. 

Although the random walk in the SWLA accurately reproduces the exact conductance 
distributions, it is not in full agreement with the statistical properties of the different trans- 
mittances. As an example, Fig. 10 shows the behavior of (T a ) = ^ b (T& a ) versus L/£ for 
both the exact numerical results (continuous lines) and the random walk (dashed lines). 
The disagreement could be associated to the use of an approximate expression (Eq. (4.17)) 
for the mean free paths. For small lengths compared to the mean free paths, the average 
reflectance is given by (see also Eq. (3.34)) 

(Rat) = + ■ ■ ■ • (4-30) 

tab 

We could then have obtained the different mean free paths £ a b for all modes by performing 
a linear fitting of the numerical results to Eq. (4.30). However, as long as the energy is not 
very close to the onset of new propagating channels, we found that the numerical mfp's are 
well described by Eqs. (4.16) and (4.17) within the numerical accuracy. The discrepancy 
could then be associated to the limitations of the SWLA. The generalization of the random 
walk method beyond the SWLA is in progress. 

In summary, we have implemented a numerical method to obtain the statistical properties 
of the transport coefficients using the diffusion equation derived in this work. We have 
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10 20 30 

L/l 

FIG. 10: Surface disordered waveguides: Average transmittances (T a ) as a function of Ljt. 
The agreement between results based on the numerical solution of the Schrodinger equation (mi- 
croscopic calculation; continuous lines) and the random walk simulation of the diffusion equation 
in the SWLA (dashed lines) is not as good as for the conductance distributions. 

extensive numerical evidence of the suitability of our model to describe the statistics of wave 
transport in disordered waveguides. It is worth noticing that our model exactly reproduces 
the conductance distributions obtained from the microscopic model even though this one 
contains as many evanescent modes as needed to perform the calculation in an exact manner. 
The only parameters needed to obtain the statistics of any transport coefficient are the mean 
free paths £ a b, as it is implied by the diffusion equation, all the statistical properties being 
fixed at any length once all £ a b parameters are fixed. 

V. CONCLUSIONS AND DISCUSSION 

The central result of the present paper is the Fokker-Planck equation, Eq. (2.10), which 
describes the evolution with the length L of a disordered waveguide of transport properties 
which can be expressed in terms of the transfer matrix M of the system. 

Our starting point is a potential model in which the scattering units consist of thin 
potential slices (taken as delta slices for convenience) perpendicular to the longitudinal 
direction of the waveguide, the variation of the potential in the transverse direction being 
arbitrary. A statistical law for the potential slices is specified, as detailed in Sec. Ill B 1: 
in particular, the parameters of a given slice are taken to be statistically independent from 
those of any other slice, so that we are dealing here with the situation of uncorrelated (at 
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least in the longitudinal direction) disorder. Our result is obtained in the so called dense- 
weak-scattering limit, denoted by DWSL in the text, in which each potential slice is very 
weak and the linear density of slices is very large, so that the resulting mean free paths 
(mfp's) are fixed (see Eq. (3.38)). The statistical properties of a building block (denoted 
by BB) of length SL, say, are first derived; the BB is then added to a waveguide of length 
L to obtain a composition law, from which the diffusion equation is eventually derived. In 
the DWSL, the statistical properties of the BB, and hence of the full system, depend only 
on the mfp's which, in turn, depend only on the second moments of the individual delta- 
potential strengths. Cumulants of the potential higher than the second are irrelevant in the 
limit, signalling the existence of a generalized central- limit theorem (CLT): once the mfp's 
are specified, the limiting equation (2.10) is universal, i.e., independent of other details of 
the microscopic statistics. 

One important characteristic of the present analysis, compared with previous ones, is 
that the energy of the incident particle is fully taken into account, a consequence being that 
the generalized diffusion coefficients appearing in the diffusion equation (2.10) depend on 
the wavenumber k of the incident wave and on the length L. 

The diffusion equation (2.10) for expectation values is very difficult to solve, the main 
reason being explained in the text, right below Eq. (4.6). The original DMPK equation 
[15] for the probability distribution of certain parameters of the transfer matrix was solved 
exactly for the unitary symmetry class only [38], whereas for the evolution of expectation 
values arising from that same equation for a large number of open channels, N ^> 1, an 
iterative procedure was developed to find the result as an expansion in powers of 1/N [10]. 
In the present case, in Sec. IV A we have been able to solve Eq. (2.10) exactly for N = 1, but 
only for a few particular observables: the solution is in excellent agreement with the results 
of a microscopic calculation. However, not even for N ^> 1 have we been able to develop an 
analytic iterative procedure like the one we mentioned above; even numerically we have not 
succeeded in developing a method to solve Eq. (2.10). We have thus tackled the problem of 
extracting information from the analysis of the present paper from a different point of view, 
based on the study of the BB itself, which was shown to have universal statistical properties. 
First, we should remark that the BB is useful not only as an intermediate step to obtain the 
diffusion equation; it is interesting as a physical system in itself, i.e., a slab. In the paper we 
obtained its statistical properties up to order 5L only, with some extension to order (SL) 2 . 
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In principle, although it represents a tedious task, the procedure could be carried on to at 
least a few more powers of SL. A similar expansion was performed in an earlier publication 
[20] . Second, the BB was used in Sec. IV B to develop the method that we called "random 
walk in the transfer matrix space", which was essential for the numerical analysis based 
on the results of the present work. The results reported in that section showed excellent 
agreement with the corresponding microscopic calculations. Efforts towards an analytical 
and/or numerical treatment of the diffusion equation (2.10) itself would be very important. 

In Sec. HID we develop the short- wavelength approximation, denoted as SWLA in the 
text, which bears resemblance to the geometrical optics limit studied in optics. The results 
of this approximation allow making a connection with some of our previous work, in which 
the energy did not appear explicitly in the analysis. We should remark that the numerical 
results of the random walk in the transfer matrix space reported in Sec. IV B were performed 
within this approximation. 

In the analysis presented in this paper, the presence of evanescent modes for a single slice 
appears in the effective potential (u r )ab that occurs in Eq. (3.16) and is used to construct 
the open-channel transfer matrix; the effective potential takes into account transitions to 
evanescent modes. Our statistical law is thus postulated for the matrix elements of the 
effective potential. However, as we mentioned in Sec. II around Eq. (2.2) and in Sec. IV B 3, 
the transfer matrix for a sequence of scatterers was constructed multiplying open-channel 
transfer matrices, i.e., ignoring the presence of evanescent modes in the combination law. 
Nonetheless, the final agreement with microscopic calculations is very good. An important 
question for future investigation is thus to understand the effect of evanescent modes when 
combining subsystems to form the whole waveguide. 

In the potential model developed here the property of time-reversal invariance is satisfied 
and the treatment is also restricted to scalar waves. In the language of random-matrix theory, 
we are dealing with the orthogonal symmetry class, or /? = 1. For possible applications 
to electronic systems, it would be interesting to extend the analysis to the unitary and 
symplectic cases, (3 = 2 and (3 = 4, respectively. 

As explained in the Introduction, in earlier publications (like Refs. [10, 15]) the notion of 
maximum entropy in conjunction with a number of physical constraints played an important 
role in selecting the distribution of the BB: in a way, that selection captured the features 
arising from a CLT. We think that it would be very interesting to investigate the question 
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whether the results presented here can be obtained within such a framework. 

Finally, since the results of our model have been compared successfully only with micro- 
scopic computer simulations, we think that it would be very challenging to measure these 
same quantities in the laboratory, in order to make comparisons with real-life experiments. 
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APPENDIX A: SOME PROPERTIES OF THE TRANSFER MATRIX 



1. Transfer and Scattering Matrices 



The transfer matrix M is closely related to a perhaps more familiar object: the scattering 
or S matrix, that relates incoming to outgoing waves: 



(Al) 




The iV-dimensional blocks of the M matrix (Eq. (2.1)), 



M 





M 12 




a j3 


M 21 


M 22 




7 5 



(A2) 



are related to the reflection and transmission matrices r, t for left incidence and r', t' for 
right incidence as 



r = -r^, t' = 5' 1 
t = (a*)" 1 , r' = /36~ 1 . 
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(A3a) 
(A3b) 



The physical property of flux conservation (FC) requires the S matrix to be unitary (SS^ = 
1) while the M matrix must satisfy the pseudounitarity condition 



M f S z M = £ z . (A4) 

This is the only condition that M satisfies in the unitary, or (3 = 2, case. If, in addition, 
the system is time-reversal invariant (TRI), i.e., in the orthogonal case (3 = 1, we have the 
extra condition 

M* = S X .MS X , (A5) 
where XI 2 and S x have the structure of Pauli matrices: 

1 \ fo l\ 

, = . (A6) 

\° _1 / V 1 °/ 

Eq. (A5) implies 

M 22 =(M n )*, M 21 = (M 12 )% (A7) 

so that in Eq. (2.1) only the two blocks M 11 and M 12 , or a and /?, need be considered. The 
relation with the reflection and transmission matrices in the TRI case is now 

r = -(a*)" 1 /?*, t' = (a*)" 1 (A8a) 
t = (at)" 1 , r' = {3{a*)-\ (A8b) 

In this TRI case, the matrix e for the BB, defined in Eq. (2.4), must satisfy the relations 

e* = S X £S X , (A9) 



so that 



e 21 = (e 12 )* (AlOa) 
5 22 = (e 11 )*. (AlOb) 



Introducing the notation 



1 = 2, 2 = 1, (All) 
the TRI relations (A10) can be written as 

it = (4)*- (A12) 
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2. The transfer matrix in terms of independent parameters 



FC and TRI symmetry imply that the actual number of independent parameters of a 
transfer matrix is 2N 2 + N. This fact can be taken explicity into account by writing M in 
a polar representation as [10] 



M 



u 

u* 



VTTx V\ 
Vx VTTx 



v 

v* 



(A13) 



where u, v are arbitrary N x N unitary matrices (each contributing N 2 parameters) and A 
is a real diagonal matrix with N arbitrary non-negative elements. 
Another useful representation was introduced by Pereyra [27] 



M 



e ih o 
e~ ih * 



V 



a/1 + 7/?7* 

rj* y/1 + rfr] 



(A14) 



where h is an arbitrary N x N Hermitian matrix (thus contributing N 2 parameters) and rj 
is an arbitrary N x N complex symmetric matrix (with N 2 + N parameters) . 



3. Multiplicativity property 

Suppose we start with a system having a transfer matrix Mi and enlarge it by adding, 
on its right hand side, another system having a transfer matrix the transfer matrix of 
the combined system is simply given by 

M = M 2 M X . (A15) 

The role of closed channels (evanescent modes) is briefly discussed in the text, around 
Eq. (2.2). This multiplicativity property is extremely useful and is extensively used in 
the analytic study along the present work. However, from a numerical point of view, the 
successive multiplication of TVf-matrices is unstable because the pseudounitary group, to 
which the transfer matrices belong, is not compact [10]. This property leads to numerical 
instabilities as the norm of the transfer matrix elements can grow without limit. In contrast, 
numerical methods based on the scattering matrix can be very effective and stable [24]: if 
we have the scattering matrices S\ and S2, the scattering matrix S corresponding to the 
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composition of the subsystems (1) on the left and (2) on the right is given by 



r = n + ti x r 2 Qt x (A16a) 

t = t 2 Qh (A16b) 

t' = (t) T (A16c) 

r' = t 2 Qr[t' 2 + r 2 (A16d) 

where 

Q = (I N - rir 3 ) _1 (A17) 
and denotes the N x N unit matrix. 



APPENDIX B: EVANESCENT MODES AND THE EFFECTIVE POTENTIAL 

In this appendix we define the effective potential for a delta slice that was introduced in 
Eqs. (3.14) and (3.16). 

Consider a problem admitting N open channels and N' closed ones. We shall eventually 
be interested in the limit N' — ► oo. The total number of channels will be denoted by 
Nt = N + N'. It will be convenient to define projection operators P and Q (with P + Q = I) 
unto open and closed channels, respectively, i.e., 

N 

P = Yl l*a> <Xa\ , (Bla) 

0=1 
N T 

Q = Yl \Xa)(Xa\, (Bib) 
a=N+l 

where \xa) represents the "transverse" state defined in Eq. (3.3). The most general solution 
of the Schrodinger equation on either side of the scattering system contains: 

i) incoming- and outgoing- wave amplitudes for all the open channels. We denote by 

~ (2) 

a p the A^-component vectors of incoming open-channel amplitudes on the left and right 
of the system, respectively, while bp^ , bp 2 ^ denote the corresponding outgoing open-channel 
amplitudes. 

ii) "outgoing" closed-channel amplitudes, denoted by the A^'-component vectors Vq , b^ 1 , 
on the left and right of the system, respectively: these are the components that decrease 
exponentially at infinity. The A^'-component vectors 5,q , Oq represent the "incoming" 
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closed-channel amplitudes, i.e., the components that increase exponentially at infinity. In 
order to have a normalizable (in the Dirac delta-function sense) wave function, closed chan- 
nels can only give an exponentially vanishing contribution at infinity, so that the components 
dg\ dg\ which we keep for convenience in the following equation, will eventually be set equal 
to zero. We shall also use the notation otpp = PaP, etc. The wave amplitudes on the two 
sides are then related by the "extended transfer matrix" [10] as follows: 



7(2) 
Op 




®pp a PQ 


fipp Ppq 


T(2) 

°Q 




OlQP OLQQ 


Pqp Pqq 


af 




1PP 1PQ 


Spp 5pq 


a {2) 

a Q 




1QP 1QQ 


$QP &QQ 



if 



>p 



(B2) 



Here we are using the notation of Ref. [10], which was developed in terms of incoming- and 
outgoing- wave amplitudes, both for the 5" matrix and for the M matrix. This results in 
an asymmetry in the notation in the two vectors appearing in Eq. (B2). Perhaps a more 
common notation expresses the M matrix in terms of waves that travel to the right and to 
the left, giving a more symmetric definition. 

The extended transfer matrix of Eq. (B2), which will be denoted by M, contains four 
Nt x Nt matrix blocks. When we set, as we already mentioned, the amplitudes o!q = Oq = 
and consider, as given data, the 2N amplitudes cip\ bp , , we obtain a set of 2(N + N') 
equations in the same number of unknowns: ap 2 ', bp*\ b^ b^ . 

The "open-channel transfer matrix" of Eq. (2.1), that relates the open-channel ampli- 
tudes on the two sides as 



(B3) 



can be obtained from the extended transfer matrix of Eq. (B2) by eliminating the closed 
channel amplitudes b^ bg\ to obtain the four N x N blocks 

1 _ 



Op 




a (3 






(2) 
dp 




7 5 




Op 



a 

7 
5 



Otpp - 

Ppp- 
Ipp — 8pq 
Spp — 5pq 



1 



Sc 



iQPi 
Sqp, 



$QQ 
1 



7qp, 



(B4a) 
(B4b) 
(B4c) 
(B4d) 
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For a delta slice centered at x = and described by the potential of Eq. (3.7), one finds 
the following extended transfer matrix: 

T AT + J_ J_ 7/ _J_ 1 1 „. 1 





a 


P 




M = 








j3* 


a* 





1 1 m 1 T 

'217k u 7k 1n i 



±^-u^- 



(B5) 



In this equation, denotes the A^-dimensional unit matrix; u is the Nt x ^Vr matrix 
constructed from the matrix elements u a b of Eq. (3.7b); K is the diagonal Nt x iVr matrix 
K a b = ^a^afe; f° r open channels (a — 1, ■ • • , iV), k a is defined as the real positive square root 
of the r.h.s. of Eq. (3.12); for closed channels (a = N + 1, ■ • ■ , Nt), we define k a = iK a , n a 
being real and positive. It will be convenient to write 



Kpp 



kp 

IKQ. 



(B6a) 
(B6b) 



Substituting the extended transfer matrix of Eq. (B5) into Eq. (B4) we find, for the blocks 
of the open-channel transfer matrix 



a 







in 



1 1 1 

UPP—FF 



li^kp 
11. 1 



-.Up P - 



(B7a) 
(B7b) 



2i \fkp ykp 

where upp is the "effective potential" referred to in the text, Eq. (3.16), and defined as: 

1 1 1 



Upp = Upp - U P Q 



N 1 



-UQQ- 



uq P . 



(B8) 



For the present scattering problem associated with the potential of Eq. (3.7), the reflec- 
tion and transmission amplitudes to open channels are to be obtained from Eq. (A8), where 
a and (3 are given in Eq. (B7) in terms of the effective potential of Eq. (B8). 

As an exercise to understand the importance of including the effect of evanescent modes 
in the construction of the open-channel transfer matrix for a given scatterer, we shall apply 
the above formulation to verify that a 2D delta potential produces no scattering at all (for 
the 3D version of this statement see Ref. [37], p. 274, Problem 3-7). 

Suppose that the potential U(x,y) of Eq. (3.7) is given by (notice that u , which we 
assume positive, is adimensional) 



U(x,y) = u 5(x)5(y - y„). 
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(B9) 



Then u(y) = uq5(u — yo) and the matrix elements u a b are now given by 

Uab = U Xa(Vo) Xb(yo) = U a U b , (BIO) 

where Xaiy) are the transverse states defined in Eq. (3.3). Notice that u a b has the property 
of being "separable" in the channel indices a and b. Defining \u) as the column vector 
[mi,m 2 , • • -ujy, ■ ■ ■un t ] T , we can write the matrix u = {u a b} as 

u=\u)(u\. (Bll) 

The projections 

\up) = P\u) = [ui, u 2 , ■ ■ ■ u N ] T (B12a) 

\u Q )=Q\u) = [u N+ i, ■ ■ ■ u N+N >] (B12b) 

allow writing upp, etc., needed in Eq. (B8), as 

(B13a) 
(B13b) 
(B13c) 
(B13d) 

We substitute the relations (B13) in the expression (B8) for the effective potential and 
obtain 



Upp = 


\up)(u P \, 


UPQ = 


\ U p)( U q\i 


uqp = 


\ U q)( U p\i 


uqq = 


\uq){u q \. 



u = \up) 

where we have defined 



1 

1 ~ to t — n — w — i \ v Qi 

In> + \vq)(vq\ 



(u P \, (B14) 



M = -^=\uq). (B15) 

The inverse appearing in Eq. (B14) can be trivially calculated due to the separable nature 
of \vq)(vq\, with the result 

1 -I N .- , (B16) 



In> + \vq)(v q \ 1 + (v q \v q ) 

so that 



\up){u P \ 

" = TT(^>- (B17) 
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The square of the norm of the vector \vq) is calculated from the definitions (B15), (B12b) 
and (BIO), and the expression (3.3) for the transverse states to find 

a =N+i yyw' K 

As an example, suppose that our 2D delta potential is located midway between the two 
boundaries of the waveguide, i.e., yo = W/2. We then find 

N+N' 1 

(vqIvq) = ^ V ; (B19) 

o=JV+l,odd V y\V> ^ 

The tail of this sum, from A 1, say, to iV + N' ^> A is 

UQ n ±^' W u N + N > 
M»Q)~W L ™~27 ln — — ■ (B20) 

a=A, odd 

Thus, as A 7 "' — > oo, (vq\vq) has a logarithmic divergence, and u of Eq. (B17) vanishes in the 
limit. As a result, the 2D delta potential inside our waveguide does not scatter at all. 

In conclusion, should we replace, in Eq. (B7), the effective potential u by the "bare" 
potential u, we would find that the 2D delta potential produces non-zero scattering. The 
presence of an infinite number of evanescent modes is essential to demonstrate that our 2D 
delta potential produces no scattering. 



APPENDIX C: THE FOURTH-ORDER TERM IN THE SECOND-MOMENT 
EXPANSION, EQ. (3.31c) 

We go back to the expression for the second moments of e for the BB, Eq. (3.31). The 
terms in the line (3.31b) vanish, being third order in the individual [e r ] l ^ b and hence in the 
potentials (v r )ab [see Eq. (3.20)]. The fourth-order terms are given in the line (3.31c): only 
the first of these three terms survives; in the other two there is no way to pair the scatterer 
indices so as to get a non- vanishing result (remember that the various e r 's are statistically 
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independent and have zero average). For the nonvanishing term we have [see Eq. (3.14)] 



ihl 
led 



(Cla) 



r>s 
t>u 



Vr)aa'{V s )a'b{ 



(VtU'(Vuh' d ) E \Vr)%V.&t)%{*u)&] (Clb) 



EE 

r>s a iRi 
t>u 

EE(^)-'(^)^)(^)^^)^} E ( Clc ) 

r>s o'/3' AV 

^^(aa',^') ^\a%(3'd) 

a'0' 



d 



d 



d 2 



(Cld) 



We now take the DWSL and find [see Eq. (3.33)] 



ton (\e {2) f h \e^] 

1WQ \ L J aft L J 



<5L 



E C(aa',cP') C(a'b, (3'd) ^ A i\',\'j,hiJ,',n'i ft. u(xr\\ 
a ,p VLai{k)£ c/3 ,(k) ^£a'b{k)£p' d {k) ^ 

(C2) 



in analogy with Eq. (3.45). We have defined 



^i^d^n{5L))= [ ! ^{x^W^W^WdxM, (C3a) 

J JK(8L) 

(C3b) 



\i+h+\'+ l x' 



TZ(SL) 



Ki 8L 



K x +K 2 
2 



— e 2 



K 1SL sm 2 



2 



Here, TZ(5L) denotes the region of integration {x > x'}, i.e., half a square of size 5L. Eqs. 
(C3a) and (C3b) are analogous to the earlier definitions in Eqs. (3.40) and (3.41). Eq. 
(C3b) is valid for K\ ^ and K 2 ^ 0. The other possibilities are 



/ \i+h+\'+n' 
1 ' 2 

/ Y+h+X'+n' \ 



/ y+h+X'+n' _ 



ATi = K 2 = 



sm 



K2 
2 



.k iS l sin^ 

2 



(C3c) 

K x = 0, K 2 ^ 
(C3d) 
#x ^ 0, K 2 = . 

(C3e) 
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We have defined 



K 2 



A aa' ',c/3' 

v-y 3,1^1 



(C4a) 
(C4b) 



We see from Eq. (C3b)-(C3e), or directly from the integral definition (C3a), that in an 
expansion in powers of SL, the leading term is quadratic in SL, i.e., 

,{SLf 



+ 



(C5) 



so that the leading term, in a similar expansion, of the fourth- order contribution to the 
second moments of e in the DWSL, Eq. (C2), behaves as 



DWS 



lim (\s^r ' [ £ (2) 1 

w<! \ L J aft L J 



SL 



E 



C{aa',cf3') C(a'b,f5'd) 



m 2 



\i+h+\'+ /jt' 



0(5LY 



(C6) 



This is the result mentioned at the end of Sec. Ill B 2 b, right above Eq. (3.48). 

In the above analysis, the BB lies in the interval (—5L/2,8L/2). If it is shifted to the 
interval (L, L + SL), Eq. (C2) is modified as 



DWS 



hm ([e^Je®] 

nu/Q \ L J aft L J 



L,5L 



X ^ C(aa', eg) Cjct'b, p'd) 

^ ^/£aa'{k)£ c ^ ^t a . h {k)t^M 



+^ 2 )(L+f ) 



(C7) 



in analogy with Eq. (3.49), while Eq. (C6) becomes 



DWS 



lim (le^fjs^] 

nwQ \ L J aft L J 



L,<5L 



E 



C(aa',cp') C(a'b,P'd) 



J^, x/La>(k)£ C f3>(k) x/l a > b (k)£f3>d(k) 

{5Lf 



i(K 1 +K 2 )L_ 



+ 0(SLf 



(C8) 



APPENDIX D: ANALYSIS OF THE GENERAL TERM OCCURRING IN THE 
CALCULATION OF THE p-TH MOMENT OF e FOR THE BB 



Eq. (3.31) which gives the expansion of a second moment of e in terms of the e^'s, 
and hence to various orders in the individual e r 's, can be generalized to an arbitrary p-th 
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moment as 

m 

Mi/" .Mp 

In the analysis that follows i/ie BB will be centered at the origin. The term under the 

summation sign in this last equation is of order ^ + h fj, p in the individual [e r ]„ 6 's, and 

hence of the same order in the potential matrix elements (v r )ab s 'i it survives only if it is of 
even order in these quantities, i.e., if n\ + • • • + fi p = 2q, say. 

Using, for convenience, a simplified notation for the indices, we express the term of order 
2q under the summation sign in Eq. (Dl) as 

( n- ■ ■ [e**c) tt = e (k • ■ ■ <n ■ ■ ■ k ■ ■ ■ «J7) 

ri>—>r n 
ti>--->t(ip 

= Yl ((^i)aai ■■•(^r M1 )a Ml _ 1 6---(^ 1 )e 71 ---(% p ) 7Mp _i/) 

ri,— ,r^ x ai,— ,a Ml -i 

7i>— >7m p -i 
Ai,— ,A M1 _i 

xh(n - r 2 ) • • ■ h{r^-x - r Ml ) 



xh{h - t 2 ) ■ ■ ■ hit^ - t„ p ) 

_ iAi,-,A M1 _ 1 i;-;wu/i,-,i' M p-in 

/ / / / oai,- ,Q Ml -ib ; en,- ,7m p -i/ ' ^ ' 

Q!i,-,a M1 _i Ai,— ,A M1 _i 

7i,->7p P -l i/i,---,y Mp _i 

We have introduced the step function h(r — s) (= 1 for r > s and = for r < s) to implement 
the correct range of summation of the scatterer indices. The function F defined in the last 
line can be read off from the equation itself; it has /ii + ■ ■ ■ + fi p = 2q pairs of upper and 
lower indices and has the structure 



•^ttl&l,- ,a 2q b 2q ~ \(^i)ii6i '■■ (^r 2q )a 2q b 2q 



,r 2q 



*fi$i;:*%Z q (ru---,r 2q ), (D3) 



where the function fa^'.'-XXS^ ' ' ' ' r2< ?) is S iven b y : 

filjl- ,i2q32q f \ _ _ To Uljl r.Q ]*2qhq 

Jai 61 , -,02,625V 1 ' '' 2 ?/ ~~ rnJai&i r^g^^q 



xTT' 2ff /^(ri-ri+O, (D4) 
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where the prime in the product sign means i ^ 1^1,^1 + ^2, ■ ■ ■ ,fii + ■ ■ ■ + /J> p -i- Two 
particular examples of the structure (D3) were already encountered earlier, in Eqs. (3.32b) 
and (Clb) above. 

The expectation value appearing in Eq. (D3) can be written in terms of the cumulants of 
the various blocks of (v r )ab s into which one can partition the product (ivjai&i ' ' ' (vr 2q )a 2q b2q- 
We first give a few examples, and then the general expression. 

i) q — 1. One can write [see Eq. (3.21b)] 

\(#ri)ai&i (^2)0262) = ^ ( a A ) a 2 b 2 )5 rir2 , (D5) 

where a second cumulant coincides with the corresponding second moment, i.e., 

k\ v) (axb 1: a 2 b 2 ) = y> 2 ) {a 1 b 1 , a 2 b 2 ), (D6) 

due to the vanishing of the first moments [Eq. (3.21a)]. 

ii) For q = 2 we have: 

(£ri)ai&i(£r 2 )a2&2(£r 3 )a3&3 (^4)0464) = [ K 2 (ai&lj ^h)^ (a 3 b 3 , a^b^Sr 

+K 2 \aib 1 , a 3 b 3 )K 2 v \a 2 b 2 , a 4 b 4 )b~ rir . A b~ r2rA + ^(ai&i, a 4 &4)«2 (°2&2, a3&3MnrA a r 3 ] 

(D7a) 

+K4 (aiftx, a 2 &2, «3^3, a 4 6 4 )^ (D7b) 

where a fourth cumulant is defined in the usual way, i.e., 

«4 ( a i^ a 2&2, a 3 6 3 , a 4 6 4 ) = ^ {aih, a 2 b 2 , a 3 b 3 , a 4 b A ) 

-^(aifei, a 2 b 2 )/j, 2 v \a 3 b 3j a 4 b 4 ) - /^(ai&i, a 3 b 3 )jj, 2 v \a 2 b 2l a 4 b 4 ) 

-/^(Mi, a A b A )iJ^ ) {a 2 b 2l a 3 b 3 ). (D8) 

In Eq. (D7), the two lines ending in (D7a) contain all possible pazr contractions, i.e., 3!! = 3 
terms altogether: this partition of 4 elements can be represented by the Young diagram: 

□□ 

(D9) 

□□ 

The last line, (D7b), contains the only possible quartet (the Kronecker delta with more 
than two indices is defined to be nonzero only when all the indices are equal). It can be 
represented by the Young diagram: 

□□□□ (D10) 
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iii) For q = 3 we have: 



^ri)oi6i (^r 2 )a 2 b 2 (^r 3 )a 3 b 3 (£V 4 )a 4 b 4 (^r 5 )a 5 b 5 (^V 6 )a G b 6 

= [^(aifti, a 2 b 2 )K^\a 3 b 3l 0464)^(0565, a 6 6 6 )5 r . ir . 2 5 r . 3r4 5 r5r6 

+ all possible combinations] (Dlla) 

+ [^4 (0161, 0262, 0363, 0464)^2 (0565, 0666)5 r . ir . 2r3r . 4 5 r6r . 6 

+ all possible combinations] (Dllb) 
+«6 (0161, a 2 6 2 , a 3 fe 3 , a 4 fe 4 , a 5 6 5 , a 6 6 6 )5 rir2r . 3r4r . 6r . 6 . (Dllc) 

The two lines ending in (Dlla) of this last equation contain all possible pair contractions, 
i.e., 5!! = 15 terms altogether: this partition of 6 elements can be represented by the Young 
diagram: 

□□ 

□□ (D12) 
□□ 

The two lines ending in (Dllb) contain all possible combinations of one quartet plus one-pair 
contraction, i.e., ( M =15 terms altogether: this partition of 6 elements can be represented 
by the Young diagram: 

□□□□ 

(D13) 

□□ 

The last line, (Dllc), contains the only possible sextet. It can be represented by the Young 
diagram: 

□□□□□□ (D14) 
It seems plausible that the particular examples given above can be generalized to arbitrary 
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q, so that we can write F^J 3 ^ of Eq. (D3) as 



aibi,— ,a 2q b2 q 



h ■ • ■ {v 



r2qj0,2qt 



biq / fa! 



UJl— ^2qJ2q 



bi,— ,a 2q b2 q 



(D15a) 



ri,— ,r2 q 



• K 2 (fl 2 g-l& 2 g— 1, a 2qb 2q )S rir2 ' ' ' ^r2q-ir 2q 



ri,— ,r 2q 



+ 



+ all possible combinations (D15b) 

K^\a 1 b u a 2 b 2 , a 3 b 3 , a^)^ (ash: cbebe) • ■ ■ 4 ( a 2g-iS-i> «2g&2 9 ) 
x5 rir . 2r3r4 <5 r5r6 • • • 5 T2q _ 1T2q + ai l possible combinations 

(D15c) 



+^ q \a 1 b 1 ,a 2 b 2 --- ^^gU..,]/^^,-,^ • (D15d) 

Again, the partition of 2g elements contained inside each square bracket can be represented 
by a Young diagram. The first square bracket ending in (D15b) contains all possible pair 
contractions, of which there are (2q — 1)!! altogether. 

Eq. (D15) shows that we can write i^'...'^* (omitting, for simplicity, the lower and upper 
indices, as well as the index (v) in the cumulants) as 

K 2 (a 1 b 1 ,a 2 b 2 ) K2(a2g-i&2g-i,a 2 g62g) , , r r r \m 
2^ J \ r 2,r 2 , . . . ,r 2q ,r 2q )a 



+ all possible combinations 



r2Ti—r2q 



(D16a) 



+d 



K 4 {aib!, a 2 b 2 , a 3 6 3 , a 4 h) n 2 (a 5 b 5 , a 6 b 6 ) K 2 {a 2q -xb 2q -i, a 2q b 2q ) 



d 



d 



x 



f (^4, r4, r4, r4, rg, r 6 , . . . ,r 2q ,r 2q ) d q 1 + all possible combinations 

(D16b) 



r±r§---r2q 



+ • 



^^ihl^M J2 f (r 2q , • • • , r 2q ) d . 



d" 



(D16c) 



>'2q 



The cumulants n 2t appearing in Eq. (D16) are defined, for t = 1,2, in Eqs. (D6) and (D8) 
respectively. 
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At this point we take the DWSL defined by Eqs. (3.38). The various fractions K^jd 1 
appearing in Eq. (D16) are finite because of the scaling assumed in Eq. (3.39). Also, the 
various summations in Eq. (D16) tend to finite integrals, and all the terms with factors 
of d "left over", i.e., from (D16b) up to (D16c), vanish. As a consequence, the cumulants 
K 4, ' ' ' K 2q, do not contribute in the DWSL: this is the central- limit theorem (CLT) that was 
discussed at the end of Sec. II, at the end of Sec. Ill B and in Sec. Ill C. The second 
cumulants k 2 enter through the various mfp's, as we see from Eq. (3.36). 

In the DWSL we thus write Eq. (D3) as 



DWS 

= C[aibl ' a *f ' ' ^'V^'aj&lS&fc 12. 34, • • ■ , 2, - 1 2,) 

Y '0161*0262 'o2 9 -l&2g-l'o2 9 &2g 

+ all possible combinations. (D17) 
We have used Eq. (3.33) and we have defined 

■ ■ f^Yr-tSl ( x 2> x 2 , ■ ■ ■ , x 2q , x 2q )dx 2 ■ ■ ■ dx 2q 

JSL 
JlZc{SL)i 

Here, i^ b {x) is the continuous version of the function \& r Yab °^ ^' c l- (3-15). The region of 
integration TZ C (SL) q arises from the appropriate step functions, Eq. (D4), that implement 
the correct range of summation of the scatterer indices, and from the type of pair contraction. 
We have added, in a symbolic fashion, in the argument of A, the information about the 
scatterer indices that have been contracted: in the above cases, the contraction was r\ = 
r 2i r 3 — r 4 ; ' ' ' 5 r 2q-i = r 2q- This information is redundant, but has beeen added for clarity. 
Eq. (D17) [inserted in Eq. (D2)] and Eq. (D18) generalize the earlier expressions (3.45) 
and (3.40). One of the "possible combinations", i.e., the one arising from the contraction 
ri = r 3 ,r 2 = r 4 , that would be indicated symbolically as 13,24, generalizes Eqs. (C2) and 
(C3). In an expansion of the integral (D18) in powers of 8L, the leading term clearly behaves 
as 

C:;S( fc ;«~( E ) s +'-' ( D19 ) 
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Consider now the particular case of an even moment of the BB e. For this purpose we 
set p = 2t in the above analysis, starting from Eq. (Dl). The lowest-order term in the 
expansion of Eq. (Dl) corresponds to fii = ■ ■ ■ = fi2t — 1 an d thus to 2q = 2t, in the 
notation introduced right after Eq. (Dl) (i.e., this term is of order 2t in the u r 's); in the 
DWSL it is found, by setting q = t in Eqs. (D17) and (D19), that its leading term in 
an expansion in powers of SL behaves as (SLY/ ^l ai b! " • ^a 2t &2t ■ Higher-order terms in the 
expansion (Dl) for the same moment are higher order in SL. The contribution to a second 
moment obtained above, Eq. (3.48), represents, for t = 1, a particular case of this general 
result. 

For an odd moment with p — 2t + 1, the first term in the expansion of Eq. (Dl), i.e., 
the one with fii = ■ ■ ■ — ^t+i = lj vanishes, because it is of odd order in the v r 's. The 
next-order terms in the expansion (Dl) have one of the \ii = 2 and all the other /Vs equal 
to 1 [for instance: \i\ = 2, fi 2 = • ■ ■ = ^2t+i — !]■ For these terms, 2q = 2t + 2, so that from 
Eqs. (D17) and (D19) we see that these terms are of order (5L) m / ^l ai bx • • • ^24+2624+2- 

The conclusion of the last two paragraphs is not altered when we translate the BB to 
the interval (L, L + SL). We have thus proven, for the moments of e, the behavior that was 
mentioned at the end of Sec. Ill B. 



APPENDIX E: SOME USEFUL DETAILS FOR SEC. IV A 



In the one-channel case, the quantity K l J^ d of Eq. (3.42) and the diffusion coefficient 
D l ^^(k,L) of Eq. (3.51), to be used in the diffusion equation (3.55), are given by: 



D ij ' hl (k, L) 



-!)' + (■ 

, \i+h+l 



■iy +1 + (-i) h + 



JK»' hl L 



k, 



2C 



(Ela) 
(Elb) 



respectively. We have omitted the channel indices, which would take the value 1. 
We can rewrite the pair of Eqs. (4.5), after multiplying the second one by e lx ° s , as 

Id A 



2 ds 

db r 

-7T- + X bi 
OS 



A + 2b r , 



-A - 2b r , 



dh 
ds 



— Xob r 



(E2a) 
(E2b) 
(E2c) 
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where 



A{s) = 2 (aa*) - 1 
b(s) = b r (s) + ibi(s) 



(af3) s e 



(E3a) 
(E3b) 



The quantities p\, P2 and ps appearing in Eq. (4.7) are the roots of the third degree 
polynomial P (p) = p 3 + x$p — 2x\ and are given by 

(E4a) 
(E4b) 

(E4c) 



Pi = 


U + V, 




P2 = 


--{u + v) 




P3 = 


--{u + v) 


% 2 



\u — V) 



[u — V 



with 



u 



Xq 

v/3 



1 + 



x 



1/2 



1/3 



3^ 

Xq 




Xq 



1/2 



When xq 3> 1, we expand u and as 

y/3 3 



4V3 105 



xq 2x 



5Xr 



v/3 



VS 3 4v/3 105 



and the roots are given approximately by 

1 



+ 



+ 



Xr 



jXn 



Pi 



o 



[ ^2 +i ^r 



r 



P2 

P3 = P 2 

We can thus write the exact solution (4.7) as a power series in 1/xq as 
A{s) = e 2s + ^ [-(1 + 2s)e 2s + e~ s cosx s] + O (\) 

Xq \ X oJ 



1 • 2 

o r (s) = — — e smxosH — ^ 



Xq 



•'Y 



3s \ 

1 — — ) e cosxqs — e 



(A 





bi{s) = — [— e 2s + e s cosxos] ^{s — l)e s smxos + o(—j. 

Xq Xq \ x q / 



1/3 



3^ 



Xq 



(E5) 

(E6a) 
(E6b) 

(E7a) 

(E7b) 
(E7c) 

(E8a) 
(E8b) 
(E8c) 
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APPENDIX F: DETAILED STRUCTURE OF THE EXPECTATION VALUES 
{F(M)) Lk IN THE SWLA, FOR THE CASE N = 1. 

The aim of this appendix is to derive the diffusion equation in the SWLA, Eq. (3.73), for 
the one channel case (N = 1), assuming, for the expectation values (• ■ ■ ) L k , a more detailed 
structure than the one that was assumed in Eq. (3.70). We shall present such derivation 
in two ways: i) starting from the BB (as in section HID), and ii) starting from the exact 
diffusion equation (4.2). 

The justification for the structure that we shall suppose comes from Sec. IV A, in which 
we solved Eq. (4.2) for the quantities (aa*) and (a/3), giving the exact result of Eq. (4.7) 
which, for large values of k£, was then expanded in inverse powers of this quantity, giving 
Eq. (4.9). Inspired by this last result we shall suppose that for an arbitrary observable 
F(M), the expectation value (F(M)) Lk has the structure: 



where the functions (F (M)) L , and (F (M))^ L depend on the energy only through the 
mfp which, just as in Sec. IV, is taken as a fixed parameter denoted by t\\ = I. In the 
second term on the r.h.s of Eq. (Fl) we have assumed, in principle, an infinite sum over m; 
however, we shall show that this sum has in fact a finite number of terms, as it occurs in 
Eq. (4.9). Just as in Sec. IVA, we have denoted the longitudinal momentum, Eq. (3.12), 
by k\ = k. For kt ^> 1, the quantities (F (M)}^ dominate the behavior of (F(M)) Lk in 
Eq. (Fl), precisely as we supposed in the "ansatz" of Eq. (3.70). 

For the particular case of Eq. (4.9) we make the following identifications: 




(Fl) 




i(l + e 2 ^), 



0, Vm, 




1 

-e 
4 



(aa*)^\ = (aa*) 



(2) 

-4,L — 







(F2) 
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and 

(«/?>i 0) = 0, 

= y~ L/l , (<*®% = 0, (afi)V L = - V^, 

(«/3>S = = 

= I (s - 3§) e-^, = = 

M?i = = -ie"^ . (F3) 

Similarly, for the expectation values that appear on r.h.s of Eqs. (4.2) and (3.58) for one 
channel, we suppose, in the SWLA, the structure 



, . OO OO v v j I 

f hl ^) L k = ((• • • f*»y L , ke) n /m ' L e^ kL . (F4) 

n=l m=— oo ^ 1 



I i hl\ \^ 

[We have used the same abbreviation as in Eq. (3.61).] Again, the functions ^(- • • ) %3 ' M y 
and /(• • .)^ m A depend on k only thorugh the mfp and, for k£ ^> 1, the quantities 

\ / m,L 

(■ ■ ■ yi h dominate the behavior of ((■ ■ .yi< hlXfl \ j n Eq_ (F4). 

I L \ I L,k 

1. Derivation of Eq. (3.73) for N = 1 Starting from the BB 

Introducing Eqs. (Fl), (F4) and the first term of (3.59) as given by (3.60) into Eq. (3.58) 
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(for N = 1), we obtain 



,imk(L+5L) 



n=l m=— oo 

oo oo 



An) 



ijhlXfi 
(K=0) 



n=l m=—oo 
\i+h+l 



2i 



(kl) n 

)1°' + E E 



. (o) 00 00 U" " V / 



hi\A {n) 



n=l m=— oo 



m >£ c imkL 



5L 



+ E 

ijhlX/i 



i+h+ i jk ( l+ sl) sin §5L 



oo oo 



+ E E 

n=l m=— oo 



m,L AmkL 



/ / 



+ 



(F5) 



Here, is an abbreviation for K l i> , as in the text. The contribution of Eq. (3.62) is not 
indicated, since it either contains inverse powers of k or is proportional to {SL) 2 . Moreover, 
we have used the expression (see Eq. (3.41)) 



A ij ' hl (k, 5L) 



. i+h+1 sin fSL 



K 
2 



(F6) 



From Eq. (Ela) we see that the possible values of K are 0, ±2k and ±4/c: thus, if K ^ 0, 
=> Ki oc kt. 

We now take to the l.h.s of Eq. (F5) all the quantities that do not have inverse powers 
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of k£, and take the rest to the r.h.s, to obtain 



(F{M))T +5L -<-<™^ 



(F{M))f- E 



\) ( \ijhl\uV ' 



(K=0) 



OO OO j 

n=l m=— oo v ' 

oo oo / 

EE E 

ijhlXfj, n=l m=—oo 
(K=0) 



(F {M))t] L+ 8L z imk{L+5L) - (F ( M ))m] L e imkL 



2C 



(kiy 



+ ^yi+h+l e ix(L+^) Sil1 f 5L 



Ki 



(...) 



oo oo 

+E E 

n=l m=— oo 



+ 



ijhlX^\^ 



m,£ imkL 



vijhl\fi\ 

I 7 



(F7) 



For W > 1 we neglect the r.h.s, so that we obtain 
(F(M))™ 4I -<F(M))f 



E 

ijhlX/i 
(K=0) 



2C 



(0) 



5L + (6L) 2 « 0. (F8) 



This last equation is the analogous of (3.71) for the case N = 1; if we make the same 
expansion as in Eq. (3.72), we finally obtain Eq. (3.73) for N — 1. 
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2. Derivation of Eq. (3.73) for N = 1 Starting from the Exact Diffusion Equation. 



We introduce Eqs. (Fl) and (F4) into Eq. (4.2) to obtain 

oo oo 

E E 



d(F(M))f , ^ ^ 1 d{F{M)t] L mkL 



dL 

+ E 

m=— oo 

E - 

ijhlXfi 
(K=0) 



E 

ijhlX/i 



n=l m=—oo 



(key 



dL 



ime imkL 



<F(M)>:v£ (F(M)>;i 



n=2 



71-1 



2f 



\ijhl\/i 



)!' + E E 



oo oo 



n=l m=— oo 



(H) n 



2^ 



\ij,hl,Xfi 



(0) 



oo oo 

E E 

71=1 m=— oo 



ijhl\{i\ 



/, 



(h£) n 



• (F9) 



If we now take to the l.h.s all the quantities that do not have inverse powers of hi, and all 
the rest to the r.h.s, we obtain 



dL 



E 

ijhlXfi 
(K=0) 



2(: 



i+h+l 



AKL 



m=— oo 
m^O 



oo oo 



ijhlXfi 



EE E 

n=l m=— oo ijhlXiM 

(K^0,K=0) 



(") 



2f 

yjhlX{i^ n 



(0) 



\ijhlXfi\ 

I L 



2f 



(kiy 



oo oo 



-E E 

71=1 772= — OO 



e 77nfcL (d(F{M)) { 2\ 



dL 



- + ^(F(M))^ n+1) 



' m,L 



(F10) 
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For W > 1 we neglect the r.h.s and obtain: 



<9(F(M))j 0) 
dL 



£ 

ijhlX/j, 
(K=0) 



i+h+l 



2C 



yjhl,\ny°) 



(Flla) 



+ 



(1) irnkL 
m,L e 



~~ etKL 



m=—oo 



ijhlXfj, 



2C 



\ij,hl,\/A 



(0) 



(Fllb) 



This last equation has two kinds of terms: the square bracket in (Flla), which is inde- 
pendent of k, and the square bracket in (Fllb), which depends on k through the complex 
exponentials. Eq. (Fll) can only be consistent if each square bracket vanishes. As a result, 
since, for N = 1, K takes on the values 0, ±2/c, ±4/c, the sum over m appearing in the second 
bracket must have a finite number of terms: m = ±2, ±4. For the soluble example of Eq. 
(4.9) we have verified explicitly that (Flla), (Fllb) and, in particular, the contribution of 
order l/k£ of the last two lines of Eq. (F10), are exactly zero. 

Starting from Eq. (4.2) we have thus obtained the diffusion equation for one channel in 
the SWLA, i.e., 



d{F (M)) 



(0) 



£ 

ijhlXfi 
(K=0) 



(") 



i+h+l 



2C 



V m hl I Mi x M l » d 2F ( M ) \ 
u \ dM iX dM h » / 

ijhlXfi 
K=0 



(0) 



(F12) 

which is Eq. (3.73) for N = 1. 

The diffusion coefficients D 1 ^ appearing in Eq. (F12) can obtained from Eq. (3.68) as: 



£,11,11 

f-jll,22 

^12,21 
^11,12 



1 

1 

Ye 
1 

2? 

511,21 



(F13a) 

(F13b) 

(F13c) 
(F13d) 
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APPENDIX G: A CENTRAL-LIMIT THEOREM FOR THE COMPOSITION OF 
BUILDING BLOCKS 



In this Appendix we obtain the statistical properties of the transfer matrix for the 
whole waveguide of length L. The waveguide is divided into P BB's of length 5L each, so 
that L = P5L. Since we can write the transfer matrix for the p-th BB as M p — I + e p , the 
total transfer matrix Mt is given by 

M T = MpMp^ ■ ■ ■ M 2 M 1 = I + 5M T , (Gla) 
p p p 

SM T = J> p+ £ e Pl e P2 +-=J2 SMp\ (Gib) 

p=l Pl>P2 A*=l 

in analogy with Eq. (3.26). Here, pi denotes the position x Pi of the center of the Pi-th BB. 
Suppose that we propose a model in which the BB's are statistically independent and that 
the various moments of the e's of the p-th BB satisfy 

= 0, (G2a) 

oo 

(feOjft (at • • • feOSfc) = E (*p) > (G2b) 

oo 

(^)lt ■ • ■ (^iss^) = E • ( G2c ) 

We shall make contact below with the physical model discussed in Sees. Ill and IV. We 
allow the functions f^ 9 ^ (x p ) and f^q+i^ (x p ) to depend on the positions x p . It is easy 
to show that in the mathematical limit P — > oo, 5L — > 0, so that = L, the statistical 
distribution of the matrix 5Mt only depends on the functions f% (x p ) appearing in the 
second moment ^(sp)^ ^pfa^b^)- Therefore, the functions f^ +r> (x p ) appearing in the 

second moment {[Sp) 1 ^ (^p)^) with a ^ and the functions fiq^ f^i+i^ ( x p) 
appearing in moments higher than the second do not play any role in the final statistics of 
8Mt- This result can be interpreted as a CLT for the composition of BB's. 

The function (x p ) was considered above as an arbitrary function of x p . Now we 
identify it with the factor that is proportional to the first power of 5L in the model discussed 
in the present paper: 

i) For the exact description, Eq. (3.50) gives 

tf\x P )=2D%$* a (k,x p ), (G3) 
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ii) For the SWLA, Eq. (3.69) gives 



if fe)=2^'3 2 . (G4) 
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